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ABSTRACT 

The problem of asymptotic expansions of Green functions in perturbative QFT is 
studied for the class of Euclidean asymptotic regimes. Phenomenological applications 
are analyzed to obtain a meaningful mathematical formulation of the problem. It is 
shown that the problem reduces to studying asymptotic expansions of products of a 
class of singular functions in the sense of the distribution theory. Existence, uniqueness 
and explicit expressions for such expansions [As-operation for products of singular 
functions) in dimensionally regularized form are obtained using the so-called extension 
principle. 
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1 Introduction 

. . . the perturbative version of the Wilson short distance 
expansion is completely clear since more than ten years. 

... A new theory is not needed. 

A referee for NUCLEAR PHYSICS B. 



1.1 The problem of asymptotic expansions in perturbative 
QFT 

Approximations and asymptotic methods pervade applications of mathematics in nat- 
ural sciences: nothing simplifies a problem more than reduction of the number of its 
independent parameters — and the idea of asymptotic expansion is one of the two most 
useful in this respect. 

In applied quantum field theory, one deals with amplitudes — or, more generally, 
with Green functions of elementary or composite operators, which depend on momen- 
tum and mass parameters. Detailed investigation of dynamics of physical processes is 
rarely possible with more than just a few independent parameters. Therefore, one nor- 
mally considers asymptotic regimes in which almost all the momenta and masses are 
much larger (or smaller) than the chosen few. Since the most informative dynamical 
framework is currently provided by the perturbation theory. Green functions are repre- 
sented as infinite sums over a hierarchy of Feynman diagrams, and the initial problem 
falls into two parts. 

First, one has to expand individual Feynman diagrams with respect to external 
parameters — masses and momenta. This is the analytical part of the problem. Em- 
pirical evidence indicated long ago that diagrams can be expanded in powers and 
logarithms of masses and momenta. For a wide class of the asymptotic regimes a for- 
mal proof of this fact was given by D. Slavnov The reasoning of was based 
on a modification of the techniques used in the proof of the Weinberg theorem (which 
established the power-and-log nature of leading asymptotics in high-energy regimes; 
for a detailed discussion of the theorem see 0). That result was extended to other 
asymptotic regimes using various techniques 0-0. 

The main difficulty with asymptotic expansions of multiloop diagrams is that for- 
mal Taylor expansions of the integrand result in non-integrable singularities. This 
indicates that the integrals depend on the expansion parameter non-analytically. The 
papers followed, with variations, the old idea of splitting integration space — 

whether in momentum or parametric representations — into regions in such a way as 
to allow one to extract the non-analytic (usually logarithmic) contribution from each 
region, eventually, by explicit integration. However, complexity of multiloop diagrams 
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exacerbated by ultraviolet renormalization made obtaining convenient explicit expres- 
sions for coefficients of such expansions unfeasible within a framework of elementary 
integral calculus. Nevertheless, expansions that involve simplest functions of the ex- 
pansion parameter (e.g. powers and logarithms) are exactly what is needed in the final 
respect in applications, and even such limited information on the analytical form of 
expansions can be useful.^ 

The second part of the perturbative expansion problem is combinatorial. The ob- 
servation that asymptotic expansions of non-perturbatively defined objects should have 
a non-perturbative form is due to K.Wilson p who also discovered that such operator 
product expansions (OPE) can be highly useful in phenomenology. W.Zimmermann 
0, using what became known as the Bogoliubov-Parasiuk-Hepp-Zimmermann (BPHZ) 
techniques, was ffist to show that it is indeed possible to construct an OPE for a class 
of short-distance asymptotic regimes. [] 

The main achievements of were, ffist, the demonstration of how the expansion in 
a "global" OPE form is combinatorially restored from terms corresponding to separate 
Feynman diagrams; second, the required smallness of the remainder of the expansion 
was proved in presence of UV renormalization. However, unlike the expansions ob- 
tained in 1^, the coefficient functions of the OPE of were not pure powers and 
logarithms of the expansion parameter (i.e. short distance or large momentum trans- 
fer) but also contained a non-trivial dependence on masses of the particles (apparently 
because masses were needed as regulators for infrared divergences). Therefore, al- 
though the results of |^ are ffimly established theorems, they fell short of providing an 
adequate basis for applications (for a more detailed discussion see subsect. below) : 
first, the expressions for coefficient functions were unmanageable from purely calcula- 
tional point of viewQ; second, infinite expansions could not be obtained in models with 
massless particles, e.g. QCD; third, the logarithms of masses contained in coefficient 



functions are — as became clear later |11| — non-perturbative within QCD in the sense 
that they cannot be reliably evaluated within perturbation theory using asymptotic 
freedom. All the above drawbacks have the same origin: a lack of the property that is 



now known as perfect factorization (we will discuss this point later on in subsect. p.2|) 



Lastly, although the BPHZ techniques may be fine as an instrument of verification 
of the results discovered by other methods, its heuristic potential turned out to be in- 
adequate: it proved to be of little help in finding new results. Indeed, both the formula 



^ cf. an early attempt to develop calculational algorithms for coefficient functions of operator 
product expansions @]. 

^ It was also realized that the short-distance OPE is closely related to other expansion problems e.g. 
the problem of decoupling of he avy particles and low-energy effective Lagrangians — see a discussion 



and references below in subsect. |2.3| ; in fact, short-distance OPE and related problems constitute a 
subclass of Euclidean regimes studied in the present paper. 

This explains why the first large-scale calculations of OPE beyond tree level |10| were performed 



by "brute force": the coefficient functions of an OPE were found by straightforward calculation of 
asymptotics of the relevant non-expanded amplitudes and then by explicit verification of the fact that 
the asymptotics have the form which agrees with the OPE ansatz. The more sophisticated methods 
of [0 were discovered outside the BPHZ framework — using the ideas described in the present 
paper. 
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of the Bogoliubov /2-operation and the OPE were discovered using heuristics that 
are foreign to the BPHZ method, while the attempts to obtain full-fledged 0P13-like re- 
sults for a wider class of asymptotic regimes (non-Euclidean, or Minkowskian regimes) 
have so far largely failed. 

A fully satisfactory solution of the expansion problem — which it is the aim of the 
present and companion publications to describe — has been found only for the class 
of the so-called Euclidean asymptotic regimes (for precise definitions see sect. in a 
series of publications [|I2|-|^, including efficient calculational formulae for coefficient 
functions of OPE. The derivation in [^-[^ employed a new mathematical techniques 
based on a novel concept of asymptotic expansions in the sense of distributions. Al- 
gorithmically, the key technical instrument here is the so-called asymptotic operation 
(As-operation). As has become clear the techniques based on the As-operation 



offers a comprehensive and powerful alternative to what is known as the BPHZ theory. 



1.2 The lesson of the i?-operation 

Bogoliubov's 1952 derivation of the fundamental formula for subtraction of UV diver- 
gences (the i?-operation) [^] provides a fine example of a highly non-trivial reasoning 
which led to a highly non-trivial result — a reasoning that has nothing to do with how 
the i?-operation is treated (i.e. "proved") within the framework of the BPHZ method. 
By examining Bogoliubov's reasoning one can exhibit the central dilemma of the the- 
ory of multiloop diagrams, and the lessons learnt thereby have a direct bearing on the 
theory of asymptotic expansions as well. 

The dilemma is as follows. On the one hand, Feynman diagrams are objects whose 
complexity increases infinitely with the order of perturbation theory. But that com- 
plexity is not amorphous, it has a structure: Feynman diagrams can be generated, 
e.g., by iterating Schwinger- Dyson equations or via some other equivalent and orderly 
procedure (see e.g. the construction of perturbation theory in from microcausality 
condition). To put it shortly, multiloop diagrams are organised in a recursive fashion. 
For definiteness, here is an example of a causality condition of the kind that was used 
by Bogoliubov: 

T[C{x)C{y)C{z)C{u)] = T[C{x)C{y)] x T[C{z)C{u)], for x°, y^ > z^, (1.1) 

which expresses the fact that chronological product of four Lagrangians is expressed as a 
simple product of T-products involving lesser number of Lagrangians, taken in a special 
order. (Each Lagrangian will correspond, in the final respect, to a vertex in a diagram.) 
On the other hand, the elements that participate in the recursion p..l| , are singular — 
their formal nature can be best described by qualifying them as distributions — and 
their products must not be treated formally: products of distributions do not, generally 
speaking, exist (if one insists that they do, one runs into UV divergences ||24 



Nevertheless, Bogoliubov's finding consists in that if all T-products of lower order 
are known, the totality of such relations defines T[C{x) . . . C{u)] everywhere except for 
the point x = y = z = u. 
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The last step in the definition of T-product can be best described using the language 
of distribution theory as a procedure of "extension of a functional" (see below sect. |]). 
Without entering into detail, we only note that at a practical level, such an extension 
consists in adding to the r.h.s. of |1.1| a counterterm localized at the point x = y = z = u. 
It can only be a linear combination of derivatives of 5{x — y)6{y — z)6{z — u). The 
number of derivatives to be added is determined by the leading singularity of the r.h.s. 
at the point x = y = z = u, which can be determined by, essentially, power counting, 
after which there are simple ways to fix the coefficients in order to ensure finiteness of 
the resulting T-product. 

The /^-operation now emerges as a straightforward iteration of the same elementary 
step in situations with increasing number of Lagrangians on the l.h.s.: addition of a 
counterterm corresponding to the singularity at an isolated point. 

The net effect of such a reasoning is, in a profound sense, "organizational" : it allows 
one to use the inherent recursive structures in order to reduce the reasoning to just one 
simple step. It is no secret (see e.g. the proof of the Weinberg theorem 0) that the 
problems involving multiloop diagrams are reduced, essentially, to a very cumbersome 
(if done in a straightforward way) power counting. A proper "organization" of the 
problem requires to do power counting only for a simple subproblem (isolated point in 
the above example) so that the solution becomes rather obvious. 

However, Bogoliubov in the 50s was not familiar with the techniques of the modern 
version of the distribution theory p3| , which explains the decision not to try to formal- 



ize the underlying heuristic derivation (such a formalization can be found e.g. in ^0 



but rather to treat the formula for the /^-operation as given and simply "rigorously 
prove" it using the simplest available method — reducing the renormalized diagram to 
an absolutely convergent integral by resolving all recursions in the framework of a 
parametric representation. Thus the BPHZ method was born. 



1.3 74s-operation: distributions and "perfect" expansions 



It turns out [|l^] that the problem of expansions of multiloop diagrams can be considered 
in a way similar to Bogoliubov's treatment of UV divergences. Since the main body 
of the present paper is devoted to explanation and clarification of this fact, only a few 
general remarks are offered here. 

First, the entire set of Feynman diagrams is structured in a recursive fashion, even 
if in more difficult problems (e.g. expansions) it may not be easy to notice the recursion 
and its relevance to the problem at hand. 

If the second part of the expansion problem (restoration of expansions in a global 
OPE-like form) is to be successful, then the expansions of individual diagrams should 
be done so as to preserve the recursive structure. The problem gets complicated if 
properties like gauge invariance are involved which connect sets of diagrams. Thus, 
an efficient handling of recursions is a key to successful organization of any work with 
multiloop diagrams. 
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A scrutiny of the problem of singularities encountered in formal expansions of in- 
tegrands from the point of view of the expansion problem in its entirety — i.e. taking 
into account the recursion aspect — reveals that the underlying fundamental mathe- 
matical problem is that of expansion of products of singular functions in the sense 
of distributions (sect. Its solution involves counterterms (similarly to the theory 
of i?-operation described above); each counterterm corresponds, as we will see, to a 
subgraph, and the underlying recursive structure of multiloop diagrams is naturally re- 
flected in the expansions, allowing simple exponentiation at the second, combinatorial 
stage |TH . 



On the other hand, in order to save effort by avoiding proving useless theorems, one 
must aim at obtaining expansions that run in powers and logarithms of the expansion 
parameter (since we know the analytical structure of expansions; otherwise one would 
have to determine it in the process). 

Now we come to the most important point: once the analytical structure of expan- 
sions is determined, a conclusion follows: if such an expansion exists, it is unique (see 
subsect. |5.2|) . An immediate corrolary is that all structural properties of the original 
collection of diagrams prior to expansion are inherited by the expansion in an orderly 
manner. This fact drastically simplifies study of gauge properties of the resulting 
expansions. 

Thus, the two key ideas of our method are: expansions in the sense of distributions 
(which allows one to make a full use of the recursive structures of the perturbation 
series) and the obligatory requirement of "perfection" of expansions to be obtained. 

One can say that, in the final respect, the two ideas offer a constructive heuristic 
framework to bridge the gap between the lowest level of the problem (the underlying 
power counting, division of integration space into subregions etc.) and its highest 
level (global non-perturbative structure of the OPE-like expansions): reasoning using 
the language of distributions and "relentlessly pursuing perfection" in expansions at 
each step, we are guaranteed that the resulting formulae will possess all the desirable 
properties. 

It is worth stressing that within the framework of our techniques the problem is 
represented as an iteration of the same simple step (construction of counterterms for a 
singularity localized at an isolated point), and one only has to ensure "perfection" of the 
expansion at this step. The property of uniqueness will take care of "perfection" of the 
complete expansion. To appreciate this it is sufficient to recall that the BPHZ theory 
(which attempts to deal with the entire problem without structuring it into simpler 
steps) failed to construct OPE in models with massless particles (which problem is 
closely related to that of constructing perfectly factorized expansions) until the solution 



was found by other methods [12 



The recipes of [13|, |14] and the multiloop calculations performed e.g. in [29| confirm 



that a good organization has its advantages. 

It remains to note that neither one of the above two ideas is specific to the case of 
Euclidean asymptotic regimes. 
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1.4 Purpose and plan 

The aim of the present paper is to explain — at a heuristic level comparable to that of the 



derivation of the i?-operation in the classic book pj] — both the current understand- 
ing of the expansion problem and the analytical ideas behind the solution |jl2[-p2[, 
presenting the results in a form that proved to be most convenient for applications in 
applied QFT. We begin by analyzing phenomenological applications in order to arrive 
at a meaningful formulation of the asymptotic expansion problem and reviewing the 
results previously known in order to compare them with ours. We identify the prob- 
lem of asymptotic expansions of Feynman integrands in the sense of distributions as a 
central one, and present its explicit solution in the form of the As-operation. In this 
paper we do not aim at achieving a complete mathematical formalization but rather 
concentrate on the key motivations, notions and ideas. We will use the dimensional 
regularization |Q (for a review see |2^) as a familiar technical framework to deal 
with various singularities one encounters in this kind of problemsj^ However, the use 
of dimensional regularization is by no means essential — a regularization-independent 



account of the proofs can be found in [20|-[22| where a fully formalized treatment of 



asymptotic expansion of products of singular functions is presented. 



In the companion publications |]T8| and [jl9| the As-operation is applied to two 
problems, respectively: deriving and studying OPE and its generalizations in the MS- 
scheme, and studying UV renormalization and related issues. In |TE| the combina- 



torics of the transition from expansions for single diagrams to OPE in a global form 
is studied. To this end the so-called As-operation for integrated diagrams is defined 
(it was first introduced in [1^). Its combinatorial structure is fully analogous to that 
of the -R-operation, which fact allows one to obtain exponentiation easily and derive 
full asymptotic expansions of perturbative Green functions in any of the so-called 
Euclidean asymptotic regimes which comprise the familiar short- distance OPE, the 
low-energy effective Lagrangians to all orders of inverse masses of heavy particles etc. 
The derivation automatically yields the convenient calculational formulae that have 
been presented earlier |TB|, and extensively used in applications (see e.g. |]2S|). It 



should be noted that the treatment of combinatorics of the ordinary /^-operation given 
r3 is, probably, one of the simplest in the literature. 



m 



In the third paper a new approach to studying UV-divergent integrals will be 
described which is based on the use of the As-operation of the present work, and a new 
representation for the /^-operation will be given which generalizes the results of sect. ^. 
(Ref. is a simplified version of the more formal text I^Tf.) It should be noted that 



one of the major difficulties in studies of asymptotic expansions of Feynman diagrams 



* The better known rigorous definition of tlie dimensional regularization is in terms of the a- 
parametric representation while our reasoning is essentially based on momentum space picture. 
Formal constructions of dimensional regularization in terms of position/momentum representations 
exist |3^, but the original definition of |Q is the most useful one in applications (which means, 
by the way, that it is this definition that should be a preferred subject of investigation; indeed, it is 
not difficult to adapt it for the purposes of formal proofs |37t , but an in-depth discussion of this point 
goes beyond the scope of the present publication). 
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used to be the presence of a subtraction procedure for eliminating UV divergences . 
The representation (see also sect. | below) trivializes the problem by reducing it to 



study of double asymptotic expansions which can be accomplished by exactly the 
same methods as in the case of one-parameter expansions. It should be stressed that 
such a simplification mechanism is not limited to Euclidean problem, it is completely 
general. An informal description of our treatment of UV divergences (presented in full 
mathematical detail in [^) will be done in [0. We will also explain how the new 



representation and the As-operation work together to produce the algorithm for cal- 
culations of renormalization counterterms and renormalization group functions known 



as the i?*-operation [30 



The present paper is organized as follows. In the introductory sect. ||, we review 
some phenomenological problems of applied QFT falling under the general heading of 
Euclidean asymptotic expansions, and analyze the properties which the corresponding 
theoretical expansions must possess in order that their applications to phenomenology 
make sense. We also review the most important earlier results on the subject and 
explain what improvements and generalizations are offered by the new theory. In 
sect. ^ we explicitly formulate the Euclidean asymptotic expansion problem which we 
wish to study. 



In sects. the basic ideas of our method are discussed. Heuristic motivations for 
studying expansions of products of singular functions in the sense of the distribution 
theory (the Master problem) are presented in sect. ^. Sect. |^ contains some basic 
background mathematical definitions. The extension principle — our key to solving the 
Master problem — is presented and discussed in sect. |[ In sect. |^ we demonstrate how 
the ideas of the extension principle work with a simple but representative example. 
Sect. H deals with expansions of one-loop Feynman integrals and introduces some ideas 
concerning the treatment of UV-divergent graphs, which are generalized in [B^ and 



19 . 



In sects. |9|-[r3|, we solve the Master problem in the general Euclidean case. First in 
sect. 1^ a system of notations is developed which allows one to compactly represent and 
conveniently manipulate Feynman integrands, without using integral representations — 
which is important if the benefits of the viewpoint of the distribution theory are to 



be used to the full. Sect. [1^ describes and classifies singularities of the formal Taylor 
expansion of singular functions with respect to a parameter, and the important notion 
of IR-subgraph is introduced. In sect. |TT| we derive the general formulae for the As- 
operation for products of singular functions, which provides the solution to the Master 
problem. Explicit expressions for the counterterms of the As-operation are derived 
in sect. |12|. Sect. ^ contains a detailed discussion of an example with a non-trivial 
pattern of singularities. 
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THE PROBLEM OF EUCLIDEAN 
ASYMPTOTIC EXPANSIONS 

2 Euclidean expansions in applications 

Since there are important aspects in which formulation of the expansion problem within 
the standard BPHZ theory had been lacking, and since in publications with emphasis 
on applications the exact mathematical nature of the problems is often not understood 
clearly, it is appropriate to review at least those phenomenological problems that are 
reduced to special cases of asymptotic expansions in Euclidean regimes. This will 
provide motivations for our formulation of the problem of Euclidean asymptotic ex- 
pansions, and give us an opportunity to review the key results on the subject as well as 
to explain the improvements which are contained in our results as compared to those 
obtained by other authors. Some important new notions will also be introduced. 



2.1 Short-distance operator-product expansions 

K. Wilson PI considered expansions of operator products at short distances of the form: 

T[A{x)B{0)] ^ Yl c^{x)Om, a; ^ 0. (2.1) 

i 

However, despite the undeniable heuristic value of the point of view of position space 
representation, it is momentum representation picture that is immediately relevant for 
analysis of phenomenological problems while the transition to coordinate representation 
somewhat obscures the Euclidean nature of the corresponding asymptotic regime. 

In fact, the generic problem which is essentially equivalent to short-distance expan- 
sions and immediately related to phenomenology, is to expand expressions like: 



J dxe'^- < T[A{x)B{Q)X{Uyi)] >o (2.2) 



as Q = -q +00. 



It should be stressed that q goes to infinity along Euclidean (space-like) directions. 
The following examples are intended to illustrate this point. 

(i) Within the QCD sum rule method (for a review see ^6[), when one studies 
the spectrum of hadrons built of light quarks, one considers expansions of the correlator 
of two local operators: 

n(g2) = i J dxe'"'' < T[A{x)B{0)] >q (2.3) 
at much greater than mf^dsy the masses of light quarks. 
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(ii) For deep- inelastic lepton-hadron scattering (for a review see |^^) the moments 
of the observable structure functions are directly expressed (see |Q) as integrals of the 
quantity 

W{Q\2pqM,pl) = I rfxe'"" < N\T[A{x)B{0)]\N > . (2.4) 



One can see from the expressions presented in |^ that the behaviour of the n-th 



moment at +oo and fixed n is determined by the behaviour of W{Q^) in the 

deeply Euclidean region 

Q^:^\2qpN\p%. (2.5) 

(This does not contradict to the popular "light-cone" philosophy where one studies 
the behaviour of the structure functions F{x,Q^), x = Q'^/2qpN, at — > +oo and 
fixed X, which is related to the light-cone expansion of the operator product. Indeed, 
reconstruction of functions from their moments is a mathematically "ill-posed" problem 



50|] and requires additional information for its solution. This agrees with the fact that 
the short- distance limit is apparently more stringent than the light-cone one. However, 
from the point of view of phenomenological usefulness the two approaches are rather 



complementary ||51| .) 

{Hi) Further examples are provided by the problem of hadronic for mf actors (usually, 
electromagnetic ones; for a general review see [Q) within the QCD sum rules approach 



Here one deals with the three-point correlator of the form (cf. Fig. 1) 

J dxdy exp i{qix -q2y) <T [J (x) J{y)j {0)]>o, (2.6) 

where the current J corresponds to the hadron, and j is the electromagnetic current. 
Denote = —qf,Q'^ = —{qi — q2Y- The case of "intermediate momentum trans- 
fer" m, 

~ Q' » (2.7) 

presents an example of a problem with more than one heavy momentum. However, this 
case still belongs to OPE problems (x and y tend to simultaneously); the complication 
due to the fact that there are more than two operators to be merged at one space-time 
point is inessential. 

{iv) The case of "low momentum transfer" 0] is described by the asymptotic 
regime 

g^»g2^rn2, . (2.8) 

This also is essentially an OPE-like problem: J{x) and J{y) are separated by a short 
distance x — ?/ — > 0. Note that here the additional spectator current j is present within 
the same T-product along with the two operators J{x) and J{y) which are to be merged 
at a point. This fact precludes a straightforward substitution of the expansion into 
|2.6| and gives rise to extra terms similar to those that appear in the example {Hi). 
However, such complications are easily treated within the formalism which we will 
describe. 
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It is convenient here to point out the difference between the Euchdean and Minkow- 



skian asymptotic regimes. Suppose the momenta in |2.6| are time-hke and one con- 
siders the so-called Sudakov asymptotic regime 

» ~ Kd... ■ (2-9) 

The crucial point is that this does not imply automatically that qi is of order m„^... com- 
ponentwise but only that there is a light-like vector pi, pf = 0, such that —pi ~ , 
while the components of pi are of order a/Q^. Such asymptotic regimes differ drastically 
from the purely Euclidean ones considered in the present paper. Although the basic 
ideas of our technique are sufficiently general to make it applicable to the Minkowskian 
case, specific implementation should take into account many new features — e.g. non- 
linearity of manifolds on which singularities of integrands are localized, etc. One can 
get a flavour of what to expect in Minkowskian case from the reviews p^, [BH]. 



2.2 Technical aspects of OPE. Perfect factorization 

Within perturbation theory, a first version of OPE was obtained by Zimmermann 0] 
in the following form: 

T[A{x + OB{x)] = J2 cnM^ m)0^,{x) + odeD- (2.10) 

i 

This result, however, despite the great theoretical significance of 0, had a somewhat 
theoretical rather than practical relation to calculations, e.g. within pQCD. (For jus- 



tice's sake, however, it should be stressed that the result |2.10| , which in any case is a 
theorem, dates back to pre-QCD times.) 

First, the above version of OPE relies on momentum subtractions for UV renormal- 
ization, while all practical calculations beyond one loop are always carried out within 
the MS-scheme introduced by 't Hooft (for a review see p6|). 



Second, for models with massless Lagrangian fields (e.g. for QCD with massless 
gluons) the dependence of the coefficient functions in |2.10| on can not be got rid of 
(this is connected with the fact that in the MOM-schemes for UV renormalization 
were used and masses in fact played the role of infrared cutoffs) — a somewhat awkward 
property from the viewpoint of practical calculations. 

Third, the applicability of the OPE to the phenomenology of QCD rests on the 
assumption that the coefficient functions are calculable within perturbation theory. 



But the coefficient functions in |2.10| depend on the masses m of elementary fields 
non-analytically and contain large log(m^^^) contributions. A renormalization group 
analysis |jTl| shows that in pQCD the renormalization-group resummation of such terms 
takes one outside the region of applicability of perturbative methods. Presence of non- 
perturbative contributions defeats attempts to extract the non-perturbative values of 
matrix elements of the composite operators by fitting the r.h.s. of the above expansion 
against experimental data — a procedure typical of phenomenological studies of deeply 
inelastic scattering, the QCD sum rules method described above etc. 
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The notion of perfect factorization — which had been lacking in, and seems to be 
inherently foreign to, the BPHZ theory — was introduced in [1^.0 Its importance has 
by now been fully realized. At the level of operators/Green functions, perfect factor- 
ization means that all the non-analytic dependence on m in the expansions like |2.10| is 
localized within the matrix elements of the operators Oi while the coefficient functions 
Ci depend on m analytically. It follows immediately that such an expansion would 
allow taking the limit m and, therefore, should be expected to be valid for models 
with massless particles. Moreover, since the masses m are no longer needed as IR 
regulators, one can expect to get rid of the dependence on N in ETTUl Furthermore, 



perfect factorization, implying that no non-analytic dependences are contained in the 
OPE coefficient functions, guarantees that extraction of the non-perturbative matrix 
elements from data is a mathematically meaningful procedure. 

At the level of individual diagrams, perfect factorization means that the asymptotic 
expansion should explicitly run in powers and logarithms of the expansion parameter. 



Such expansions are unique (see below subsect. 11.2 ) which fact is of paramount tech- 
nical importance; for example, it allows one to conclude (cf. subsects. [11. 2| , |12.5| ) that 
such asymptotic expansions commute with multiplication by polynomials and other 
algebraic operations on diagrams — the property that greatly simplifies combinatorial 
study of OPE etc. in non-scalar models; in particular, one need not worry about things 
like gauge invariance of the expansions since such properties are inherited by expan- 
sions from initial integrals in a well-defined and orderly manner Uniqueness of such 
expansions also allows one to exhibit the recursive structure of the problem of expan- 
sion of Feynman diagrams which can be effectively used in derivation of OPE (see 
below subsects. |11.3H11-^ - 



It turned out that the OPE derived directly within the MS-scheme [|T^,[] pT| , 
possess the property of perfect factorization so that: the dependence of the coefficient 
functions on masses is analytic (which was earlier shown to be a necessary condition of 
existence of the OPE in the MS-scheme 0); the expansion is a full asymptotic series; 
it is valid in models with massless particles; and due to greatly simplified dependence 
of coefficient functions on dimensional parameters, there exist algorithms for analytical 
evaluation of the relevant integrals through three loops All these properties are 



interrelated, as is stressed above and as will be seen from our derivation. 

More conventional derivations of OPE in the MS-scheme a la BPHZ were presented 
by Llewellyn Smith and de Vries and Gorishny [?D|. The latter work clarified the 



^ In the context of the theory of i?-operation perfect factorization corresponds to renormahzation 
using the so-called massless schemes among which the MS-scheme is by far the most useful one. 

^ The techniques for studying gauge properties was extended to the framework of the Euclidean 
As-operation in [ p8[ . In this respect recall that gauge invariance plays a central role in studies 
of non-Euclidean asymptotic expansions which represent a major unsolved problem with important 
applications to physics of hadrons (see ) . 

^ which result came about as a special case of the general theory of EA-expansions developed in 
||l2|-[^, Efficiency of the techniques developed therein allowed us to obtain very general 

formulae first made public in . The results of were widely discussed in the literature , |Q , 

i3. 
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connection of the results of [jT2|-[|T6| with the BPHZ method.0 



Theoretically speaking, one need not use specifically the MS-scheme to obtain OPE 
with the property of perfect factorization. Any of the so-called massless schemes (or 
generalized MS-schemes [^) will do p2 |. 



It should also be clearly understood: the most important thing about the new 
results is not the fact that, say expansions were proved "in the MS-scheme", but that 
such expansions possess the property of perfect factorization. This should not be 
surprising in view of the fact that the MS-scheme has the property of polynomiality of 



its counterterms in masses which itself is a special case of perfect factorization — 
this point is best understood in the context of the definition of generalized MS-schemes 
2^ by way of subtraction of UV-dangerous asymptotics from momentum-space 



m 



integrands, so that evaluation of those asymptotics is done using the same apparatus 
of the As-operation as used in obtaining asymptotic expansions of diagrams. This 
similarity is not so much technical as ideological. 



2.3 Heavy mass expansions 

Consider a QFT model in which all the fields are divided into two sets: light fields 
collectively denoted as ip and heavy fields $. Their Lagrangian masses are m and M, 
respectively. Assume that all M are much larger than m, and consider a Green function 
G of light fields only. It depends on M besides other — momentum, coordinate etc. — 
parameters. The problem is to study the asymptotic expansion of G in inverse powers 
of M : 

GL..,M) ~ ? (2.11) 

M^oo 

In fact the first two terms in the heavy mass expansion had been known for some 
time. The "decoupling" theorem |Q (for a review and references see [0) asserts 



that, up to 0(M~^) corrections, the effect of the virtual presence of heavy particles 
(in phenomenologically meaningful cases) is absorbed into a finite renormalization of 
the Lagrangian of the light world, whose parameters become dependent on M. This 

® Still another approach, specifically designed for achieving ultimate rigour in treatment of dimen- 
sional regularization within the techniques of the a-parametric representation is sketched in . It is 
unclear, however, how successful the latter attempt was because neither that paper nor the publica- 
tions cited therein contain explicit treatment of the analytical part of the proof (recall in this respect 
that a-parametric integrands for MS-renormalized diagrams are rather cumbersome distributions [Q 
so that advantages of using such a representation are somewhat unclear). On the other hand, the com- 
binatorial part of the proof is discussed in the original publications , |]l8j in a more straightforward 
manner due to recursive economy of the formalism of the As-operation. 

The regularization independent aspects of the original derivation have been exhaustively discussed 
in pO| , pH and p^ . Various ways to rigorously treat dimensional regularization without paramet- 
ric representations were discussed in p6| , although not much practical understanding is added 
thereby to the heuristic treatment of p3[ . Anyway, the variety of uses of dimensional regularization 
in practical calculations is such that there is little hope that everything that is being done will ever 
be "rigorously proved". In the final respect, mathematics is not about rigour. Mathematics is, first 
and foremost, about calculations. 
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fixes the form of the leading term of the M — > oo expansion. The next-to-leading 
term can be described as differing from the initial Green function by one local operator 
insertion. (This result was first obtained in [5^| by making use of Wilson's OPE.) 



More formally, let £tot be the total Lagrangian of light and heavy fields, and let < 
T(p{xi) . . . (p{xn) >o'°' be a Green function of the light fields in the world containing 
both light and heavy fields. Then the leading and next-to-leading contributions to the 
expansion have the form: 

r AO) 
< T^{xi) . . . ^{xn) >^*°* ^ z"(M) < T^{xi) . . . ^{xn) >o'" (2.12) 



+ < Tlp{xi) . ..ipiXn) 



>J-^ + o(M-^), (2.13) 



(0) 



M 

where is obtained from that part of £tot which describes the light fields, by a finite 
M-dependent renormalization of its parameters; C^J^{y) is a local operator composed 
of light fields; z{M) and g{M) are functions of M with logarithmic leading behaviour 
at M — s> oo. (The UV renormalization in all cases is implicit.) 

The notation is suggestive, and indeed, to all orders in M~^, 



where 



<T(^(a;i)...(^(x„) >^-' ~ 2;"(M) <T(^(a;i)...(^(x„) >^^«, (2.14) 

Af-+oo 



£eff = E^^""^?n(^)4?- (2.15) 



The derivation of this formula will be given within the MS-scheme, so that our results 
will contain as a by-product a derivation of the decoupling theorem and of the formula 
2.13| directly within the MS-scheme, which has been lacking. Moreover, the expansions 



we will obtain possess the property of perfect factorization similar to that discussed 
in the preceding subsection, so that the heavy and light parameters are fully factor- 
ized, i.e. z{M) and g{M) do not depend on m, but only on the heavy mass M and, 
logarithmically, on a renormalization parameter. 

The most important phenomenological applications of the heavy mass expansion 
are: the study of effects of heavy quarks on the low-energy light quark properties, espe- 
cially on the parameter Aqcd (for a review see ^^); the low-energy quark Lagrangian 
of the electro-weak interactions Il57| ; evaluation of muon contribution to the electron 



anomalous magnetic moment [£8[. The latter problem is interesting as an example of 
a problem where a very high precision — and, consequently, taking into account many 
terms in the expansion — is required. 



2.4 Effects of heavy masses on OPE 

The QCD sum rules approach to studying the heavy quark bound system [^] 
presents an example of a problem where a short- distance OPE is modified by presence 
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of a heavy mass. Here one should study the correlator p.3| but for such an asymp- 
totic regime that the squared mass of the heavy quark were of the same order of 
magnitude as Q^. (This point was stressed in |5^ .) 

In the context of deep-inelastic lepton-hadron scattering this asymptotic regime 
was studied in pO . 

It follows from our general results that the expansion in this case has the form 
similar to ordinary OPE: 



eq.2.3 



J2ciiM,Q)<0,iO) >o, 



(2.16) 



but with coefficient functions depending on the heavy mass M. To our knowledge, an 
all-order derivation of the expansion p.l6| has heretofore been lacking. 



2.5 Light-by-light scattering. Linear restrictions on heavy 
momenta 

Consider the problem of deeply inelastic light-by-light scattering with both photons 
deeply virtual Here one studies the four-point correlator of hadronic electromag- 



netic currents with the pattern of external momenta as depicted in Fig. 2a, in the 
asymptotic regime 

g?~g^»m2,,. {Q',=qf). (2.17) 

A new feature in this case is that along with the general momentum conservation an 
additional linear restriction is imposed on the heavy momenta (cf. Fig. 2b). Motivated 
by this example, we will allow arbitrary restrictions of this kind in our general treatment 
of Euclidean asymptotic expansions. 

Note that the Green function to be expanded may not exist at the restricted mo- 
menta owing to infrared singularities. On the other hand, if it does, its expansion can 
be obtained via our technique in the same way as in the non-restricted case. However, 
the Q-dependent coefficient functions of the expansion of the non-restricted Green 
function (cf. the one corresponding to Fig. 2b) will in general be singular at the re- 
stricted momenta and have no relation to the coefficient functions for the restricted 
case. In other words, restricting heavy momenta and performing expansion are in 
general non-commuting operations. 

Allowing restrictions on the heavy momenta has no bearing on our formalism for 
expansion of a single multiloop Feynman graph but only on the subsequent combina- 
torial analysis of the global structure of the resulting expansions of Green functions as 
a whole. 
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2.6 Contact terms 



In our treatment of Euclidean asymptotic expansions we will be considering Green 
functions as distributions with respect to the external moment a.0 Apart from funda- 
mental meta- level arguments (finite energy resolution of measuring devices etc.), there 
are more specific motivations for this. 

The first one is of a technical nature: we will in any case make an essential use 
of the theory of distributions in order to expand integrands of multiloop Feynman 
graphs so that considering Green functions as distributions over external momenta 
brings uniformity to our argument (cf. subsect. [4.2| below) and is therefore inherently 
natural. 



Another — rather typical — motivation comes from the QCD sum rules method [pS 



The starting point here is a spectral representation for the correlator |2.3| , e.g. 

/o "'" s + Q 



U{Q') = I ds^^^^. (2.18) 



An important element of the techniques of is the so-called Borelization procedure 



which is applied to both sides of p.l8| in order to suppress contributions from large 



s in the spectral sum on the r.h.s. of p.l8| . It is assumed in ||3^ without proof that 
Borelization procedure commutes with the asymptotic expansion of n(Q^) at 
+00. However, this problem can be avoided (at the cost of some complication of 
formulae) by the following modification of the arguments of P^. Consider the relation: 



F{s/X) = fd-Q^^ = const H d..-'^^^. (2.19) 
^ ' ' J s + Q^ Jo s + x ^ ^ 

One can find — e.g. using tables of integral transforms — pairs of f{Q'^) (related to the 
Bessel functions Ju{VQ^)) and F{s) (related to the McDonald functions K^{s)) such 
that: 

(i) eq. |2.19| is satisfied; 

[a) ^p{Q'^) is smooth in Q everywhere including Q = 0; 

{Hi) ^pIQ"^) is bounded from above by some negative power of at oo; 

(iv) F{s) decreases at s ^ oo faster than any power (namely, as exp(— y^) 
times a negative power of s). 

Then using p.l8| and |2.19| one gets an exact sum rule: 



//•oo 
d''Qip{QyX)U{Q^) = dsp{s)F{s/X). (2.20) 

The weight function F on the r.h.s. provides as good a suppression of larger s 
contributions as the simple exponent used in |^ . At the same time, the l.h.s. is the 

^ The fact that expansions in this case run in powers and logs was obtained in ||^ , without exphcit 
expressions for the contact terms. There are many other examples in the literature where distributions 
and asymptotic expansions appeared together — see e.g. eq. (46) in fzot. 
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value of the distribution Il{Q^) on the test function ip{Q'^/\) (the fact that ^{Q"^) does 
not decrease at +00 as fast as is required of the Schwartz's test functions 

is immaterial since we do not have to deal with arbitrary tempered distributions but 
with a specific one, Il{Q^), whose behaviour at large is sufficiently good to make it 
well-defined on the (p{Q^)). 

The next step in the recipe of would be to use the expansion at Q"^ +00 for 
Il{Q^). Let us rewrite the l.h.s. of p.20| as 



A^/^l d''Qy^{Q^)Il{XQ^). (2.21) 



Now it is clear that Il{Q^) should be expanded not as a function but as a distribution 
"in the sense of distributions" . Let us explain the meaning of this phrase in our specific 
case (cf. subsect. ^.1|) . 

The asymptotic expansion of n(Q^) considered as an ordinary function of at 
Q'^ — * +00 in the usual sense can be written as 

n(AQ2) ^ 5]c,(AQ2)0„ (2.22) 

where are coefficients independent of Q^, and Ci(Q^) are ordinary functions, e.g. 
Cj(Q^) = Q^*-""^- One can not integrate |2.22| against a test function because, first, 
Cj(Q^) may be too singular at Q — * 0; second, there is no guarantee that termwise 
integration of the r.h.s. results in a correct asymptotic expansion of the integral on the 
l.h.s. 

On the other hand, for an expansion 

U{\Q') ^ Y.^''i>^Q')0:>, (2.23) 

A^oo ., 
t' 

to hold in the sense of distributions, all cf/''*''(Q^) must be well-defined distributions 
over the whole range of Q, and its termwise integration against any test function tp{Q) 
must generate a correct asymptotic expansion at A ^ 00 of the A-dependent function 

/(A) = J d^gn(AQ2)^(Q). (2.24) 

The relation between 2.22| and 2.23 can be established by noticing that both expansions 
allow integration with test functions localized in small neighbourhoods of the points 
Q ^ 0. It follows that cf/^*'' at Q 7^ should be either zero, or it should coincide with 
some Ci] in the latter case O- = Oj. Therefore, 

Y.4''''iQ)o', = T.^fiQn + T.^^AQn', (2.25) 

i' i i' 

where Acj/(Q) are localized at Q = 0, which implies that they are constructed of 6{Q) 
and its derivatives, while cf{Q) are distributions which by definition coincide with 
Ci{Q) at Q 7^ 0. For example, if 

n(Q)^i^Oi + o(Q-^), (2.26) 
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then the expansion in the sense of distributions will have the form 

1 .R 



n(g):^^ — ^ Oi + 5(g)oi + o(g-^), 



where 



1 



R 



lim 



Q2 + e 



(2.27) 



(2.28) 



with chosen so as to render finite the integrals of |2.28| with arbitrary test functions. 
Note that the expression |2.27| is unique — cf. below subsect. |5]^ — so that a change of 
by a finite constant is compensated by an appropriate change of O'. And, of course, 
one can use any other regularization in p. 28 . 



We see that the expansion in the sense of distributions differs from the expansion 
in the usual sense by ^-functional contributions (compare |2.26| with |2.27| and p.28| ). 
Such contributions can be conveniently denoted as contact terms. Ref. contains 
an example which demonstrates that such contact terms can be numerically significant 
in applications. 



Our formulae for general Euclidean asymptotic expansions ||T6| 
expressions for the contact terms. 



18| include explicit 



3 The formal problem of Euclidean expansions 

Motivated by the examples considered above, let us now present a formal description 
of the general problem of Euclidean asymptotic expansions, following W^. We will 



stay within the framework of perturbative QFT ||2J] and assume that the MS-scheme 
5H] is always used for UV renormalization. 



3.1 Green functions 

There are two types of fields collectively denoted as cp (light fields) and $ (heavy fields), 
with Lagrangian masses m (some of which may be zero) and M. Let hj{xj) be local 
monomials built of (p, and Hi{yi) built of both if and (One could allow h to be 
also built of $ but such cases seem to never occur in applications.) Denote the full 
Lagrangian of the model as 

Aot(¥',*,m,M). (3.1) 

Consider the Green function 

G{Q, M, k, m,fi) = J (J[ dxj) {j\ dyi^ exp i Qiyi + ^ kjXj (3.2) 



/ 3 



X < vac\ 



RT\[{Hi{yi)\{h,{xj)\\vac >^'°\ (3.3) 
I i 
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where R is the UV i?-operation and n the renormahzation parameter; T denotes the 
chronological product. It is assumed that the correlator is evaluated within the model 
described by £tot; in the functional notation: 

< vac\T[. . .]\vac >^t°t^ J dipd^ [. . .]expi J dxCtotix). (3.4) 



The Green function 3.2 as defined should be proportional to the 5-function which 



expresses the overall momentum conservation, but we will ignore it, just assuming that 

Y.Qi + Y.kj = o- (3.5) 

I j 

The momenta Q and masses M will be called heavy while k and m, light. It will 
be assumed that the ratio of scales of light over heavy parameters (denoted as k) is 
vanishing, and G is to be expanded with respect to k. More precisely, we wish to expand 
in K the expression 

G{Q,M,Kk,Km,^) - ? (3.6) 

A remark is in order. We have implicitly assumed that the renormahzation param- 
eter n belongs to the heavy parameters. However, since we always use the MS-scheme 
for UV renormahzation, each Feynman graph contributing to the Green function G is 
a polynomial in log fi. Therefore, if we take fi to be proportional to k, the additional 
dependence on k can be easily taken into account. At the level of Green functions, the 
renormahzation group can be used to get rid of the ^-dependence of fi. The net effect 
is that it is irrelevant whether we consider as a heavy or light parameter. Having this 
in view and using the homogeneity of G with respect to its dimensional parameters: 

G{Q, M, Kk, Km, /i) = k'''"'''G{Q/k, M/k, k, m, /i/zt), (3.7) 

one sees that the problem |3.6| is equivalent to studying expansions of G in Q,M oo, 



which is more readily associated with phenomenological problems. The form |3.6| of the 
expansion problem is more convenient for our purposes. 



3.2 Perfect factorization 

The most important requirement which asymptotic expansions should satisfy in the 
framework of perturbative QFT is that the dependence on heavy and light parameters 
should be fully factorized (see the discussion in subsect. ^3). Thus, if the expansion 
has the formQ 

G{Q,M,Kk,Km) ~ y2ci{Q,M)Gi{Kk,Km), (3.8) 

t 

then not only Gi should be independent of Q and M, but also q(Q, M) of k and m. 

It has become customary to refer to Ci as to "coefficient functions" and to Gi as to "matrix 
elements" — the terminology inherited from the theory of short-distance OPE. 
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At the level of individual diagrams, the requirement of perfect factorization means 
that individual terms on the r.h.s. of should be be proportional to a power of k 
and be polynomials of log k. This automatically implies that q are sums of powers and 
logs of Q and M, while Gi, of k and m. In a more general context, one should aim at 
obtaining expansions |3.8| in such a form that the scaling properties of, say, Ci{Q,M) 
with respect to Q and M be as simple as possible.^ 

The requirement of perfect factorization is completely universal and not limited to 
the case of Euclidean regimes. 



3.3 Euclidean regimes 



The problem |3.6| is yet too general. However, the examples considered in sect. 
demonstrate that it will remain sufficiently meaningful if we limit ourselves to the 
special case of Euclidean asymptotic expansions. Namely, we assume that 



< 0, for each subset C of /. (3.9) 




The problem ^]6| with the restriction ^[9] can be called Euclidean asymptotic expansion 



problem. It allows a full and explicit solution ||12||-[p2 



The meaning of the restriction |3]^ is as follows. Perform the Wick rotation of all 
integration momenta in Feynman diagrams into the Euclidean region. Then only the 
light external momenta k will remain non-Euclidean. However, since one performs 
expansions with respect to at = 0, the fact of k being non- Euclidean is inessential 
from technical viewpoint. Therefore, from the very beginning one can assume one 
works in a purely Euclidean theory. This drastically simplifies the geometry of IR 
singularities^ and, in the final respect, allows one to obtain an explicit and complete 
solution of the general Euclidean expansion problem. 

Note also that it is for the Euclidean version of the general expansion problem that 
most clean phenomenological predictions can be made within QCD where asymptotic 
freedom allows one to perform renormalizat ion-group improved calculations of the Q- 
and M-dependent coefficient functions of OPE-like expansions. 

From now on, we are going to work in a purely Euclidean theory, so that the 



condition 3.S is satisfied automatically. 



A priori the form of dependence on the expansion parameter may not be known. In more general 
contexts — e.g. beyond the framework of perturbation theory — the language of dilatation group and 
its irreducible representations might be useful here. 

For example, consider the singularity generated by the massless propagator l/p"*. In Minkowskian 
region it is smeared over the light cone ~ while in Euclidean region it is localized at the point 
p = 0. 
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3.4 Linear restrictions on heavy momenta 



For the heavy momenta on the l.h.s. of p.2| the following restriction due to momentum 
conservation holds: 

J2Qi = -Kj2k, ^ 0, (3.10) 

, . K — >u 

I 3 

i.e. the sum of all heavy momenta is not itself heavy. A less trivial example of linear 
restrictions is provided by the light-by-light scattering problem (subsect. p.5| ). 

In general, one can allow any number of linear restrictions of the form 

Y,ciQi = 0{k) (3.11) 

to be imposed on the heavy momenta. In such a case we may assume that there exist 
linearly-independent momentum variables Qi which are independent of k, and each Qi 
is expressed as a linear combination of Qi and k with coefficients independent of k : 



Qi = Qi{Q,k). (3.12) 



The restrictions p.lO|, p.ll| should thus be automatically satisfied. 



One should take care not to impose restrictions such as would make the Green 
functions to be expanded ill-defined. However, even if the restricted Green function 
develops a singularity at k = but is otherwise well-defined for all k 7^ 0, then the 
expansion problem still makes sense (the singularity will then show up as a contribution 
proportional to a non-positive power of k) and can be treated by our methodsQ. 

Among all possible linear restrictions, the so-called natural restrictions are of im- 
mediate physical interest. Such restrictions have the form 

^ Qi = 0{k) - 0, (3.13) 

where is a subset of the set of all /, and £q, for different a are pairwise non- 
intersecting. In other words, the external lines corresponding to heavy momenta are 
organized into non-intersecting bunches, and within each such bunch the heavy mo- 
mentum conservation holds separately. 

One could also impose similar linear restrictions besides conservation on the light 
momenta k. We will not do this, however, because in all interesting cases the final 
formulae will not be affected thereby (unlike the case of restrictions on heavy momenta). 



3.5 Contact terms 



As was discussed in subsect. 2.6, it is natural to consider Green functions as distribu- 



tions over momentum variables; the expansions in the sense of the distribution theory 

It should be noted that in this case the singular dependence on k, in the expansion will be localized 
within the matrix elements G, in |3.8| while c, will be insensitive to what value k, is set to. 
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are characterized, from the pragmatic viewpoint, by presence of contact terms propor- 
tional to 5-functions of heavy momenta Q in q on the r.h.s. of p.8| . 

More specifically, let F{Q) be a test function of the independent variables Qi. Then 
what one has in fact to expand is the expression 

^(l[dQ^)F{Q)G{Q{Q,Kk),M,Kk,Km,^^) ~ ? (3.14) 

I 

for arbitrary F. This means, first, that all Cj in |3.8| are distributions and, second, that 
|3.8| retains its asymptotic nature after termwise integration against F{Q). 

When there is no need to take into account contact terms, it is sufficient to consider 
Q as parameters analogous to heavy masses. Then we will say that one deals with the 
simplified expansion problem without contact terms. One should only be careful to fix 
Q at "non-exceptional" values at which the Q-dependent factors of the expansion are 
smooth in Q. The exact criteria here are of little practical usefulness. 

Note that in most cases of phenomenological interest the Green functions are inte- 
grable (although singular at some points) functions of the light momenta k. Therefore 
for the sake of simplicity we do not introduce test functions for k; it will be sufficient to 
assume that k are fixed at some non-exceptional values. Normally, if the Green func- 
tion and all terms of the asymptotic expansion are integrable functions of k, then the 
expansion allows termwise integrations against fc-dependent test functions; or, to put 
it formally, the operation of asymptotic expansion commutes with such integrations. 

To summarize: we work within the perturbative QFT, always use the MS-scheme 
for UV renormalization, consider arbitrary models, and wish to obtain full asymp- 
totic expansions for arbitrary Green functions for the class of Euclidean asymptotic 
regimes. We allow linear restrictions to be imposed on the heavy momenta of the 
Green functions, and consider and expand the latter in the sense of distributions. In 
the expressions to be obtained, the heavy and light parameters must be fully factorized. 



BASIC IDEAS 
4 Why expand products of singular functions? 

The central notion of the theory of As-operation is that of asymptotic expansion in the 
sense of distributions. In this section we present motivations for studying expansions of 
products of singular functions in the sense of distributions. Although we have assumed 
to work in Euclidean theories, our reasoning here remains valid in the most general 
Minkowskian case as well. 
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4.1 Local study of multiloop integrands 

We start with studying expansions of an arbitrary multiloop Feynman diagram G which 
can be taken to be one-particle irreducible. Let G{p, k) be its integrand (prior to UV 
renormalization) where p is the set of all loop momenta, and we have also explicitly 
shown the dependence on the expansion parameter k. The final goal is to obtain an 
expansion in k of the expression 



R j dpG{p,K), (4.1) 



where R is the UV i?-operation, while integration over p is performed in infinite limits. 
A scrutiny reveals, however, that the starting point should be the expansion not of the 
integral [4.1| but of the integrand considered as a distribution over p. 



Indeed, whatever prescription for the i?-operation were chosen, one would not be 
relieved of the necessity to study contributions from finite regions of the integration 
space. Therefore, it is natural first to replace the i?-operation by a cut-off at large 
integration momenta. This allows one to strip the expansion problem of the collateral 
complications due to UV renormalization — anyway, many diagrams are not UV diver- 
gent at all. The cut-off can be chosen in the form of a smooth function H{p/A) equal 
to zero at all sufficiently large p : 



J dpG{p,K)H{p/A). (4.2) 



A straightforward approach to expanding such integrals in k would be to first expand 
their integrands. However, the formal Taylor expansion of G{Q, k) in powers of k 
results, as a rule, in non-integrable singularities. For example, for a scalar propagator 
with a light mass m = 0{k,) one has: 

11.1 , , 

m'- + ... (4.3) 



p2 _j_ ^2 p2 pi 



The singularities of the formal expansions like those on the r.h.s. of |4.3| will be called 
Euclidean infrared — or simply IR — singularities. They correspond to the so-called soft 
singularities which appear in studies of IR divergences in QED and QCD, while collinear 
singularities have no analogue in the Euclidean case. 



The non-integrable singularities on the r.h.s. of |4.3| preclude substitution of 
into |4.2| . Emergence of such singularities also means that the integrals depend on 
K non-analytically. From the technical point of view, it is our major task to develop 
systematic algorithms for treating such singularities in this and more complicated cases 
and to construct correct asymptotic expansions starting from the formal series like 



The traditional approach to the problem is to split the integration region into two 
parts: a neighbourhood O of the point p = 0, and the complement of (9, and then to 
extract the non-analytic part of the expansion coming from the integral over O by e.g. 
explicit integrations or using special tricks. Although it is the most straightforward 
way to establish the power-and-log nature of the expansion (cf. flj, 0), it is practically 
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impossible to obtain explicit expressions that were suitable for restoring the OPE- 
like combinatorial structures out of pieces coming from different diagrams. There is, 
nevertheless, another approach based on the ideas and notions of modern mathematics, 
namely, the theory of distributions [^, which we now proceed to describe. 



Our key observation is that the same propagator |4.3| can appear within different 
diagrams and even in different places in the same diagram, so that it is natural to 
regard it as a distribution, the rest of G{p, k) playing the role of a test function. The 
same is also true of groups of factors in G{p, n) containing more than one propagator, 
and of G{p, k) itself. More precisely, let (9 be a small region of the integration space. 
Let G"^™^(p, k) be the product of those factors from k) whose formal expansions are 
singular within O. Denote the rest of G as G'^^^. There exist many Feynman diagrams 
leading to the same G^™^ but different C'^^'s. Therefore, G^™^ can be regarded as 
a distribution in C, while G'^^^ plays the role of a test function. (Such a reasoning 
leads to a very important localization property of the As-operation considered below 
in subsect. |11.3| .) 



One can also observe that besides troublesome propagators like [4.3| , there are 
factors — e.g. propagators with a heavy mass, 

^ (4.4) 



M2 + {p-Kky' 

or vertex polynomials — whose expansion in k is harmless. Denote the product of all 
such factors in G as G''*^^ and the rest of G as C^™^. Again, G^'^^ may vary for the 
same C"^™s^ go that the latter plays the role of a distribution and the former of a test 
function. 

Summarizing our observations, we arrive at the conclusion that the fundamental 
technical problem is to learn to expand in k the integrals like 



dpip{p) G{p, k) (4.5) 

with arbitrary test functions ip{p); this is exactly what is meant by saying that one 
has to obtain an expansion of G{p, k) in the sense of distributions (for a precise 
definition see below subsect. |5Tl| ). 



The problem of studying expansions of products of singular functions with respect 
to parameters in the sense of the distribution theory will be referred to, in view of 
its paramount importance within our formalism, as the Master problem. It plays the 
crucial role in our theory, and all mathematical difficulties of analytical nature are 
concentrated in it. 



4.2 Remarks 



{i) If if in ^]5| is localized within a small region (9, then one can study local structure 
of the expansion of G{p, k,) in each O. Then an expansion valid for all p can be easily 
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obtained using standard techniques of the distribution theory. Such a locahzation trick 
proves to be a powerful instrument in obtaining expansions of products of singular 
functions — see below sect. |ll] and . 



(ii) The transition from localized (f to (f = 1 requires an analysis of asymptotics 
of the corresponding distributions in the UV region p oo. For example, as will be 
shown below, studying UV renormalization of the graph G is equivalent to expansion 
of G{Xp, k) considered as a distribution in p 7^ 0, at A — > 00. One can show that this is 
a special case of the Master problem with a special choice of the expansion parameter. 
Studying how asymptotic properties of an expansion are affected by taking the limit 
p 00 reduces to studying double As-expansions which is done in essentially the same 
manner as in the case of expansions with respect to one parameter ||22[ . 



{in) The expansion problem as formulated above requires to consider Feynman 
diagrams as distributions in heavy external momenta. Such a formulation is very 
natural in the context of the Master problem. Indeed, if the diagram is integrated over 
the heavy momenta Q with the test function x{Q)y then instead of the integrals H?^ 
one should expand integrals of the form 

J dpdQ^{p)x{Q)G{p,Q,K). (4.6) 

Combine both sets of momenta into one: p' = {p,Q), consider G as a distribution in 
p', and expand the expression 



dp'y^'{p')G{p',^). (4.7) 
One can see that, essentially, there is nothing new here as compared with [4.5|. Choosing 



^'{P^ Q) = ^(p) xiQ)^ we return to Kb 



4.3 Analogy with UV renormalization 

As was pointed out in there exists an analogy between the problem of UV diver- 



gences in Bogoliubov's interpretation and the problem of singularities in expansions 



like |4.3| : in both cases the divergences spring up as a result of formal manipulations 
with objects which by their nature are distributions. Therefore, Bogoliubov's solution 



of the UV problem provides an insight into how the divergences in expansions |473 
should be handled. 



Let us reason as follows (for definiteness consider [4.3|) . An important observation 
is that the expansion |4.3| (with 4-dimensional p) allows (i.e. retains its meaning as an 
asymptotic expansion after) a termwise integration with any test function which is 
equal to zero in any arbitrarily small neighbourhood of the point p = 0. Therefore, 



the only way to improve upon the formal expansion is to add terms localized at 
p = to the r.h.s. So, it is to be expected that an expansion valid in the sense of 
distributions will look like 
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where c(m) is a scalar coefficient. 

The role of the new "counterterm" on the r.h.s. is two-fold. First, it should provide 
an infinite contribution to counterbalance the divergence due to the singularity of 
and ensure integrability of the r.h.s. as a whole against arbitrary test functions. 
Second, it should ensure the asymptotic character of the expansion. Note that the 
dependence of c(m) on m can not be trivial, for there is no other place for the expected 
non-analyticity in m to show up. 

Thus, heuristically speaking, our problem is to obtain explicit expressions for the 
coefficients c(m) in and in more complicated cases. 

5 Asymptotic expansions of distributions 

Although asjTiiptotic expansions and distributions did appear simultaneously in vari- 
ous contexts in the literature, it seems that no systematic study of the corresponding 
mathematical notions has ever been undertaken. Below are presented background 
mathematical definitions and some results of a general character. Since we are not 
aiming at attaining full mathematical formalization in the present paper, some tech- 
nical details are omitted and others treated at a heuristic level. A substantial formal 
treatment of the subject can be found in [pO[] . 

5.1 Definition 

Let F{p, k) be a distribution over p depending on the parameter k. We say that the 
series 

~ (5.1) 

represents an asymptotic expansion in the sense of distributions if a termwise integra- 
tion of with any test function ^{p) results in a correct asymptotic expansion in the 
usual sense, i.e. 

+ o(«:^). (5.2) 

This implies that each is a well-defined distribution over p. We say in such cases 
that the expansion |5.1| allows integration with arbitrary test functions. 

We will usually have to deal with expansions such that each is a polynomial in 
log K : 

In 

Fn{p,K) =Y.hg\ Fn,i. (5.3) 
This form of K-dependence will be called soft (see also below subsect. |5.3|) . 



dpF{p,K)ip{p) ~ X! dpFr,ip,K)ip{p) 
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5.2 Uniqueness 

If an expansion in powers and logarithms (like |5.1| , |5.3|) exists, it is unique. Indeed, 



integrating [ST^ with test functions we get asymptotic expansions in the usual sense. 



for which the uniqueness is an elementary fact ||6^. It is important to understand 



that there is no point in discussing uniqueness unless the form of dependence on the 
expansion parameter is explicitly fixed, as in |5.1|, [57 



The property of uniqueness is extremely important. First, it allows one to immedi- 
ately obtain useful technical results like the localization property of the As-operation 
for products of singular functions (subsect. |11.3 ) or its commutativity with multiplica- 



tions by polynomials (see below subsect. |12.5|) . Second, it implies that the expansions 
which we will derive cannot, in a sense, be improved upon. 



5.3 Regularization 

It is often convenient to introduce a regularization into both sides of |5.1|. We assume 



that it is controlled by a parameter e, and taking it off corresponds to e ^ 0. Denote 
the regularized version of the distribution / as {fY, so that 

fip) = lira {f{p)y (5.4) 



on each test function ip. Then instead of |52, one can write 

{F{p,^)r Y.'^'^i^nip,^)}'. (5.5) 



We stress that the expression pT5| need not be a true asymptotic expansion for e 7^ 0, but 
must only become such upon taking the limit e ^ 0. Still, the dimensionally regularized 
expansions that we obtain for the standard Feynman integrands are asymptotic ones 
even for e 7^ 0. This peculiar fact, however, seems to have no practical implications. 



Within the dimensional regularization, p.3| should be modified. Indeed, the depen- 
dence on K takes the form 

{F^ip,K)y = Y.'^''{KJ% (5.6) 

i=0 

which will also be qualified as soft. 



An important fact is that although the expansion ^]6| as a whole does have a well- 
defined finite limit at e ^ by construction, the quantities {F^j(p)}% in general, do 
not: if one integrates such an object with a test function then the resulting expression 
will contain pole singularities as e 0. It is, of course, always possible to rearrange 
the r.h.s. in such a way as to make each term have a finite limit at e — >■ 0. To this end 
it is sufficient to form special linear combinations of the braced expressions. However, 
such a representation is of little practical interest, at least as long as one stays within 
the calculational framework of dimensional regularization and perturbative QFT. This 
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is because in applications one has to perform calculations not with individual diagrams 
but with their sums. Summations over full sets of relevant diagrams result in con- 
siderable cancellations between them, especially in gauge models like QCD, so that 
the efforts spent on separate treatment of individual diagrams will be lost. On the 
other hand, there exist formulae ||13[, connecting the phenomenological quantities 



that one needs to calculate in the final respect, with the corresponding full sums of 
diagrams. Such formulae are compact and, as experience shows I^Qf, very convenient, 
which explains our use of expansions in the dimensionally regularized form. 



5.4 Av-dependent test functions 

As is clear from the motivations presented in subsect. |4.1| , "natural" test functions that 
emerge after localization of the expansion problem, depend on the expansion parameter 
K, SO that a typical problem is to expand expressions like 

dpip{p,K) F{p,k). (5.7) 

Strictly speaking, knowing how to expand F{p, k) in k in the sense of distributions 
allows one only to expand integrals with cp independent of k (or with trivial — e.g. 



polynomial — dependence). Nevertheless, under certain conditions the expansion of 5.7 



can be obtained by substituting the expansion ^]T| for F and a similar expansion for Lp 
and reordering the terms in the product in ascending powers of k. There are two such 
conditions 

[i) The first one is of little practical significance and is included for completeness' 
sake. Both F{p, k) and all Fn{p, k) must be distributions in the precise meaning of the 
word, i.e. (a) be linear functionals defined on all Schwartz test functions and (b) possess 
certain properties of continuity. However, there is no practical chance to encounter a 
functional satisfying (a) but not (b), as such functionals can only be proved to exist, 
by making use of the notorious Axiom of Choice, outside the realm of practical math- 
ematics. There would be no need to mention (b), but there has not yet been enough 
time for the alternative Axiom of Determination due to Myczelsky and Steinhaus (see 
e.g. p3) to purge the mathematical textbooks of pathological counterexamples. 



{ii) The really important condition imposes restrictions on the asymptotic expan- 
sion of the test function. Namely, if (p is expanded as 

i^{p, k) ~ ^ (/?„(p, k) (5.8) 



Our discussion at this point is simplified in order to avoid somewhat cumbersome technicahties — 
the inequahties describing continuity properties of the distributions involved, the corresponding semi- 
norms etc. Complete details can be found in feol. Sufhce it to mention that the remainder term 



of expansions of expressions like 5/7 is studied using representations of the form (for simplicity we 
consider here expansions to leading order in k) LpF — ifoFo = ipa{F — Fq) + {ip — ^a)F, where sub- 
script denotes the expansion of the corresponding term to o{k'^). Analyzing the non-trivial second 
addendum on the r.h.s. one arrives at, loosely speaking, the two conditions in the main text. 
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(where the dependence of ipnip, i^) on n is soft, as in |57^) , then the expression 

[^{P^k)-Y. K-^r.{p.K)]/K'' (5.9) 
n<N 

should tend to zero as k — > in the following sense. The expression in the square 
brackets must be non-zero only within a compact region K of the integration space for 
all K 7^ 0; the expression as well as its derivatives in ]?, must tend to zero uniformly 
with respect to p. 

For the Euclidean expansion problem within the standard perturbation theory one 
has to consider only K-dependent test functions of the form (cf. [7.43| ) 

(p{p,K) = tpip) f{p,K), (5.10) 

where ip{p) is a Schwartz test function, while f{p, k) is a rational function of p and k 
and such that all its singularities are localized at points p where il){p) = 0, for all n. It 
is not difficult to understand that for such ip{p, n) condition {ii) is satisfied. 



6 The extension principle 



Having developed a suitable language, we can turn to the expansion problem proper. 
Our purpose now is to offer a very general and abstract framework for studying asymp- 
totic expansions of distributions. We will introduce some notions and present simple 
but important propositions under a general heading of extension principle which con- 
stitutes, essentially, a context in which to work with explicit problems. To make it work 
for a particular problem, specific estimates and bounds must be established, which may 
technically be the most cumbersome part of the problem. But the importance of the 
conceptual framework of the extension principle consists in the fact that it allows one 
to guess the structure of the expansions to be obtained, shows what kind of estimates 
are required, and is, therefore, a powerful heuristic tool. 



6.1 Motivations 

Let F{p, k) be a distribution in p G P (where P is a Euclidean space) depending on the 
expansion parameter k. For simplicity assume that F{p, k) is an integrable function of 
p for all K 7^ 0. 

A typical situation which we will have to deal with can be described as follows. 
One can obtain (e.g. using the Taylor theorem) an approximation]^ 

Fip,K) ~ P^(p,/s:) + o(/€^), (6.1) 

^'"^ In this section we consider not full asymptotic expansions, but approximations to a given order 
o{k^). Transition to infinite asymptotic series can also be treated in an abstract manner but there 
is not much wisdom to be gained from it, so we postpone the discussion of this point till the next 
section where an example is analyzed in detail. 
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where F^ip, n) is an integrable function for p E O where O is an open region in P, and 
has non-integrable singularities at the points from P\0; the notationf^ o{k,^) stands 
for the terms that vanish faster than — s> 0. Suppose that the expansion lOI 



IS 



such that it is vahd pointwise for each fixed p E O, and allows termwise integration with 
arbitrary test functions localized within O. It is required to obtain an approximation 
for F{p, k) to the same order o(k^) valid in the sense of distribution theory on the 
entire P. 

Let us look at the problem from a different angle. Let L be the vector space 
of all test functions v^(p), p G P, and the set of the test functions ip localized 
within O. forms a vector subspace in L. The distribution P(p, k) is, by definition, 
a linear functional on L, which we will denote as P^- On the other hand, Fm{p,k) 
generates a linear functional P^.a? defined on L^v, and the restriction of P^ to Lat can 
be approximated as 

P«,^ + o(«:^), (6.2) 
which is to be understood so that W% should be valid in the usual sense if one eval- 



uates both sides of on arbitrary vectors ip E L^. The problem is to construct an 
approximation 

P, ~ F,^^ + o(«:^), (6.3) 

which would be valid on the entire L. Therefore, one can say that the problem consists 
in extending the approximating functional defined on the subspace onto the entire space 
in such a way as to preserve its approximation properties. Its solution is provided by 
the so-called extensions principle 



6.2 The extension principle 

Let L be a linear space and P^ a linear functional on it which depends on the expansion 
parameter k. Assume there exists a subspace Ln G L and a functional P^.a? defined on 
Ln, such that Pk.at approximates P^ to order o(k^) on L^. More precisely, for each 

<P.-P.,^,(^> =o(/€^). (6.4) 

We wish to construct a functional Pk^at defined and approximating P^ to order o{k^) 
on the entire L, and coinciding with P^^Af on L. 

The more familiar notation 0{k'^) is often used in theoretical physics to denote terms behaving as 
K^, possibly with logarithmic corrections. Note that a priori we can say nothing about the behaviour 
of the remainder of the full expansion we wish to construct except that it should vanish faster then 
the last term retained. That is why the o-notation is preferable to the O-notation in our case. 

Our extension principle belongs to the Hahn-Banach type results. The classical Hahn-Banach 
theorem (see any textbook on functional analysis) considers extension of functionals preserving the 
property of being bounded by a seminorm. In our case, the extensions should preserve the approxi- 
mation property. 
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To this end, let be a set of vectors transverse to and such that each Lp E L 
can be uniquely expanded as 

^ = ^K^a + ^, (6.5) 

where ip E L^. \t is sufficient to define F^^Af on ipa, because on ip its value is known. 
Now, the desired definition is: 

<F«,7V,</'a> = <i^K,</'a> + 0(k^), (6.6) 

where we have indicated that the definition allows an arbitrariness of order o(k^), 
which should not be surprising. 

It is not difficult to check that F^^at thus defined is unique within accuracy o(k^). 



6.3 Counterterms and consistency conditions 

The above construction can be represented as the following two-step procedure. Sup- 
pose it is possible to find an extension of tv from L^r onto the entire L, which 
does not necessarily approximate on L but has, perhaps, some other nice properties. 
Let (5/3 be a set of functionals on L such that 

<5/3,V9> = 0, for all V9 G Lat, (6.7) 

and 

<(5/3,V9q,> = 5a,p (the Kronecker symbol). (6.8) 

(Instead of |6.8| , one can simply require that 5p be a maximal linearly independent set 
of functionals satisfying |6r7| .) Then F^^]^ should be representable as 



F^,N + j:^ai'^)^c., (6.9) 

a 

where c^(k) are unknown coefficients. F^at thus defined reduces to F^at on each 



G Lat, so owing to it is sufficient to require that ^]9| approximate F^ to o{k'^) on 
each : 

<[F^-F^,N],ip^> = o{K''). (6.10) 



This consistency condition [|I2| , ||15| allows one to find the unknown coefficients in 

c^iK) = <[F. - F^A ^a> + o(k^). (6.11 



It is not difficult to check the equivalence of this construction and that of subsect. |0 . 
Note that the recipe |6.9| is more in the spirit of the distribution theory, for one can 
find a correct form of the required extension |6.9| , up to some coefficients, practically 
without reference to test functions that only make their appearance in the consistency 
conditions 16.10 
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6.4 Remark 



Assume that ipa in J-lOj depend on a parameter A without violating 6^: ip^ = '^^^ 
assume that there exists a finite hmit: 



lim < 

A— >CXD 



(6.12) 



even if y?^ itself does not converge to any vector in L. Now, if taking the hmit A ^ oo 
does not take one outside the o(k^) uncertainty, then instead of |6.11| one may take: 



lim < 

A^oo 



(6.13) 



The option to take such hmits proves to be very important — see below subsect. 173 



6.5 Regularized form of the extension principle 



There is a third way to represent the extended approximating functional -Fft,Ar, which 
will be widely used in what follows. Define a regularization of F^^n, i-e. a functional 
F^ defined on the entire L which depends on an additional regularization parameter 



e and such that F^ 



N ~^ 

lim< 



on Lat as e — > 0, i.e. 



F^^ — Fi^^N , Lp> = 0, for each Lp E L 



N- 



It is easy to check that the definition |6.9| can now be rewritten as 



where 



K,N 



lim 



< 



+ o(k 



F. - f: 



iPa> + 0{k 



N\ 



Here one can also use parameterized sequences of ipa as in |6.13 



(6.14) 

(6.15) 
(6.16) 



6.6 Summary 



The procedure of constructing an approximating functional now is as follows. One 
has the linear functional F^ on L which is to be expanded in k, and one is given 
another functional F^^n which approximates to a given accuracy on many but not 
all vectors from L. The key point is to identify the maximal subspace L^r on which F^^at 
approximates F^. After this is done, the rest consists essentially in adding to F^^at a 
linear combination of functionals which vanish on Ljv, with coefficients determined 
by the consistency conditions. 



In the context of subsect. the fact that Sa must be zero on L^v can be interpreted 
as that the corresponding distributions are localized in the complement P\0 of the 
open set O. 



Let us now turn to an example. 
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7 Example: expansion of (p^ + z^^) ^. 

The abstract considerations of the previous section provide us with a general approach 
to constructing expansions of distributions. Now we wish to show how it works in 
practice. There are also some important details — e.g. the transition from approxima- 
tions to a given order to full asymptotic expansions — which can be more conveniently 
explained using an example. 



7.1 Formal expansion 

Let P be a 4-dimensional Euclidean space; its elements will be denoted as p. Consider 
the scalar propagator 

F{p,K) = ^^. (7.1) 

We wish to obtain its expansion for k ^ valid in the sense of the distribution theory. 
First consider the formal expansion to order o{k^) : 

^ -^-^'-.+o{^'). (7.2) 



p2 + «2 p2 

It holds pointwise for any p ^ 0. Moreover, it allows integration with any test 
function which is equal to zero in some small neighbourhood of p = 0. But most 
important is that although it does not allow integrations with arbitrary test functions 
owing to the non-integrable singularity in the 0{k'^) term on the r.h.s., still it allows 
integration with test functions (f which only satisfy the condition (f{0) — 0, apparently 
because then (p{p) ~ p, p — > 0, and the singularity is effectively suppressed. 



7.2 Approximation properties of the formal expansion 

More precisely, one has to prove that 

j d'^p(p{p) 



p^ + p^ 



P 



0[K 



The simplest way to do this is as follows. One splits the integral as 



|p|<ck; 



\>CK 



(7.3) 



(7.4) 



For IpI < CK, one rcscalcs p — > rep, uses the scaling properties of the bracketed expression 
and the fact that (^(p) ~ p, p ^ 0, and finds that 



J \p\<CK 



(7.5) 
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For \p\ > CK, one uses the Taylor theorem to estimate the bracketed expression by an 
expression similar to the first discarded term in |7.2| , namely, const x n'^/p^, estimates 
vip) by const X IpI and integrates over ck, < \p\ < M < +00 where M depends on (p. 
The result is 

/ ^ 0{K'^\ogK) = o{k'^). (7.6) 



The desired result |7^ follows immediately from |7.4H7.6| . 



We have discussed the proof in detail because the same steps are performed in the 
most general well. 



Now we are going first to implement the recipe of subsect. |6.3| . After this is done, 
we will cast the results into a regularized form corresponding to subsect. using the 
dimensional regularization. 



7.3 Intermediate operation R 



Before proceeding further, let us construct a distribution R[p ^] defined on all test 
functions and such that the equality 



1 

7 



R 



1 

p^ 



(7.7) 



is valid on the test functions ip defined above. An explicit expression for such a distri- 
bution can be found, if one fixes any one test function ip{p) such that ip{0) = 1, and 
integrates both sides of |7.7| with (pip) — (p{0)ip{p). Then one finds: 



d^R 



p4 



Vip) = d^p{^ Ivip) - <^(0)^/'(p)] + c6{p)p{p) 



where 



1 d^p R 


" 1 " 







i'{p) 



(7. 



(7.9) 



is a single constant to be fixed — which can be done arbitrarily — in order that R[p~^] 
be fully defined. Introducing a regularization into the r.h.s. of |7^, one gets a represen- 
tation of R[p~^] "in terms of infinite counterterms" : 



R 



p4 



lim 



-A +zj{p) 



(7.10) 



where {■■■}" denotes the regularization, e.g. 



p4 



p^ + e^ 



1 

~4 



e(H>e)- etc.. 



(7.11; 



where we have introduced a convenient function 



<d{x) = 1, if the logical expression x is "true". 



(7.12) 
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0, otherwise. (7-13) 

It might be helpful to note that the distribution R[p^^] defined above is similar to the 
so-called "+" -distribution, 



dx 



dx. 



X 



(7.14) 



which is widely used in the QCD parton model calculations. 



7.4 Consistency condition 



According to the recipe of subsect. we should first redefine the r.h.s. of [7.2| in such 



a way as to make it well-defined on all test functions. To this end we replace p in 
by R[p~^] without changing values of the r.h.s. on and obtain an expression 



+ o(k' 



(7.15) 



which holds, as [7.2| , on (p, but unlike [7.2| , the r.h.s. is well-defined on all (p{p) 



The only distribution that can be added to the r.h.s. of |7.15| without violating 
the asymptotic character of |7.15| on is S{p). The immediate conclusion is that the 
expansion of (p^ + k^)~^ to order o(k^) in the sense of the distribution theory can only 
have the form 



1 



p4 



+ 4(k)5(p) + o{k' 



(7.16) 



and is defined up to a single numeric- valued (non-analytic) function of k,,c^{k,). 



In order to determine c^(k) 



recall the prescription of subsect. |6.3| . One notes that 
any test function <f{p) can be represented as 



ip{p) = ip{0)ip{p) + ip{p), 



(7.17) 



where ip{p) has already been defined and <^(0) = 0. On such ip 7. It degenerates into 
|7.2| and, therefore, is valid. Hence, in order that [7.1(j| be valid on all test functions, one 
only has to ensure that it holds on ipip)- Let us take this as a definition of c^(k) (a 



"consistency condition", cf. |6.10|) . Then, integrating both sides of |7.16| with ip{p), one 
gets: 



d'^pipip) 



P^ + 



pz 



+ K^R 



+ 



(7.18) 



Note that c^(k) is determined at this stage only up to o(k^). This agrees with the 
fact that ip has been chosen arbitrarily, because a change in if) preserving the condition 
■0(0) = 1 will change c^(k) by a term of order o(k^). 
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7.5 Dependence of counterterms on k 



Although eqs. 7.16 , [7.18 provide a correct expansion of the propagator in the sense of 
the distribution theory to o(k^), one further important refinement is possible in the 
spirit of subsect. |6^ . It is natural to try to find a ip{p) such that would simplify 
as much as possible. To this end, notice that the square-bracketed expression in [7. IS 
is bounded by const x K^/p^ at p — > oo. Therefore, c^(k) will be changed by a finite 
o(/t^) contribution if one simply takes ip{p) = 1. More precisely, let ip{p) = $(p/A) 
and let A oo. The resulting expression for c^(/t) is: 



(fp- 



1 



P^ + 



- — + k'R 



p4 



(7.19) 



The distribution R[p ^] has a simple scaling behaviour: 

F{\p) = X-^ [F{p) + const X log X6{p)] 



(7.20) 



(which can be easily deduced from the definition 'LS). Therefore, replacing p ^ np in 
|7.19| and using |7.2CI| , one gets: 



(7.21) 



where Ki and K2 are constants independent of k. The fact of a simple dependence of 
c^(k) as defined in |7.19| on k singles out this definition from the family |7.18| . 



7.6 The remainder of the expansion 



Another way to check the validity of the expansion |7.16| with c^(/t) given by [7.19| , is to 
notice that the remainder can be represented as: 



P^ + 



1 

p^ 



- c^i^)S{p) 



(7.22) 



d'p Ivip) - ^(0)] 



1 



1 



K 



(7.23) 



p^ + p^ p^ 

(By the way, this representation verifies that the r.h.s. of [7.16| is independent of the 
arbitrary constant in the definition of R[p~^] [7.8| .) 

To check that the r.h.s. is o(k^) one splits the integration region into two subregions: 
\p\ < Kk and \p\ > Kk, with K appropriately chosen. For the first subregion one 
represents ip{p) — ip{0) as Y.Pi'fiip) where all ipi{p) are smooth everywhere including 
p = 1^ , and then performs the scaling p — > Kp. In the second region, one bounds the 
first factor in the integrand on the r.h.s. by a constant, and uses the Taylor theorem 
to estimate the second factor. Then one easily obtains the desired estimate. (For more 
details see 



20|. 
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7.7 IR-counterterms as UV-renormalized integrals 



Let us discuss the expression \1.19[ Introduce again the cut-off $(p/A) exphcitly into 



the integrand of [7.19| . Then one can rewrite |7.19| as: 



hm 

A^oo 



/ 2 



(7.24) 



and the hmit is finite by construction. 

The first term on the r.h.s. is the vacuum average of the operator T[ip'^{x)] (without 
normal ordering), while the rest play the role of its UV-counterterms. 

We take note of this fact for two reasons. First, the coefficients of IR-counterterms 
analogous to in more complicated cases will also admit simple interpretation as 
contributions to Green functions with local operator insertions. Second, reversing 
the course of thought, one can view |7.19| as a peculiaiQ representation of an UV- 
renormalized one-loop Feynman graph in terms of subtraction of asymptotics of the 
integrand at large integration momenta. Such a representation will play an important 
role in our further study of EA-expansions, effectively reducing the case of UV-divergent 
diagrams to that of UV- convergent ones (cf. below subsect. |0 and pT|). 



7.8 Dimensionally regularized form 



The form [7.16| , |7.19| for the expansion is not the only one possible, and, definitely, not 
the most convenient one. If one uses a regularization, then, taking into account |7.19| , 
one has the following expression instead of |7.16| and [7.19|: 



1 



lim 



— K 



p4 



+ cI{k)6{p) + o{k^ 



where 



Cn K 



d^p 



P^ + 



p'i 



(7.25) 



(7.26) 



Much simpler expressions emerge within the dimensional regularization. (Note that 
in this case one regularizes also the expression to be expanded as well as those terms of 
the expansion which actually do not need to be regularized, but this can do no harm. 
For simplicity of notation we will omit the limit signs.) Thus, instead of [7.25H7.2^ one 
gets: 



where 



Co 



p2 -|- 
k) = J d^p 



— - K^— + cq{k)5{p) + o{k^ 



pz. 



p^ 



1 



p^ + 



K 



pz 



p4 



J p'^ + K 



(7.27) 
(7.28) 



or natural — depending on one's prejudices. 
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(we take e = (4 — D)/2). The terms on the r.h.s. vanished owing to the fact that the 
dimensional regularization preserves the formal scale invariance, so that e.g. 



1 



P 



,2a 



X 



D-2a 



d^'p 



1 



P 



2a 



0. 



Another consequence of this fact is: 



Co(k) = [n'^Y''^ ^Co(l) 



(7.29) 



(7.30) 



i.e. the dependence on k takes the form discussed in subsect. 5.3. 



Of the three forms of the expansion ( 7.16 and 7. IS ; 7.25 and 7.26 ; 7.27 and 7.28 ), 
the most important are: the first one since it does not use artificial regularizations 
and contains only explicitly convergent terms, and the third one which owing to its 
simplicity is most convenient in applications. Note, however, that the representation 
in terms of dimensional regularization |7.27| , |7.28| is not devoid of some mysticism: 
indeed, the IR divergence of the term on the r.h.s. of 7.27 is compensated by 



a divergent part of Cq{k) which, judging from the r.h.s. of |7.28|, has an UV origin. 



However, one should firmly remember that that is due to the fact that the dimensional 
regularization nullifies certain types of integrals (cf. |7.29|) and formally preserves scale 
invariance which establishes a connection between IR and UV contributions. Therefore, 
it becomes possible for an expression to contain both types of divergences balanced so 
as to render the whole finite. In more conventional regularizations there occur only 
compensations of divergences of the same nature (cf. |]22[| ). 



7.9 Complete expansion 



Now let us turn to higher terms in the expansion |7.16| . To derive them, instead of 
following the recipe of the extension principle, we will use a trick motivated by the 
representation of the remainder of the expansion in the form of |7.22| .P^ The advantage 



of it is that it provides a short-cut for the reader familiar with the dimensional regu- 
larization to arrive at the final formulae. However, this trick cannot be extended to 
more complicated situations of non-Euclidean asymptotic regimes. 



Consider the expression 



j d^pA{p)B{p) = o(k='^), (7.31) 



where 



1 ^ 1 



B{p) = ■p(p) - T^Mr), (7.33) 



19 



An essentially similar trick was used in [pO| 
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a;^ = 2(iV-l), (7.34) 

and 

T.^¥'(p) = E V"^(0) (7-35) 



is a partial Taylor series. The expression 7.31 will turn out to be the remainder of an 
expansion to order o(/€^^). 

The fact that the r.h.s. of |7.31| is indeed 0(^2^) can be explained as follows. B{p) 
is 0{p'^'^~^^) at p ^ and 0{p'^'^) a.t p —>■ 00, and is independent of k. A{p) is o{k'^^) 
for p 7^ 0, has a non-integrable singularity 0(p~2{^+i)) at p — > 0, and is bounded by 
Q^i^2{N+i) ip2{N+2)-^ at p — > 00. One can see that the non-integrable singularity of A{p) 
at p ^ is suppressed by B{p), while the growth of B{p) at p ^ 00 is not rapid 
enough to spoil integrability of A{p) at infinity. A more formal proof proceeds along 
the lines of subsect. fTGl 



There are two remarks worth making here. First, the interplay of asymptotics in 
7.31| is only possible owing to the uniformity of {p'^ + k'^)~^ with respect to simultaneous 



scaling in p and k. Second, the expression on the l.h.s. of [7.31| contains no cut-offs, and 
is therefore devoid of arbitrariness. 

To represent the final result in a more convenient form, introduce a full basis of 
homogeneous polynomials Va{p), such that 

VaiXp) = Al"'P„(p), (7.36) 

and the dual set of derivatives of the 5-function : 

d''pVa{p)6/3{p) = 6a,0 (7.37) 

(cf. |6.8| ). Then the Taylor expansion can be represented as: 



t^Mp) = E ^"(p) [/ dp's^ip'Mp') 



\a\<uj 



(7.38) 



For the example under consideration, one can take 

Va{p) ^ p'", Up) d^'^ip), a = 0, 1, 2 ... , (7.39) 

owing to the 0(-D)-invariance. 

Finally, introduce dimensional regularization into the (convergent) expression on 
the l.h.s. of |7.31| . Then integrate it using fr29[ Simple transformations using fr38| lead 
to the following final result, wherein / d^pip{p) has been omitted from both sides: 



p2 -|- ^2 



N 1 

EHT^+ E c„(«:)5.(p) + o(.-) 

n = " IrvKl.Iivr 



(7.40) 
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where 



1 



'11 



(7.41) 



p2 + 

and c^Tv is given by |7.34| . One can easily find from [7.41| , wliicli order in k eacli coun- 
terterm belongs to, and one can see that the additional dependence on k in each order 



in K is soft in the sense of subsect. 5.3 



To avoid misunderstanding, we stress once more that the expansion [7.40| is to be 
understood in the following sense (cf. the remarks in subsect. |5.3| ). One should first 
integrate both sides termwise against a test function, and take off the regularization 
(i.e. take the limit D —>■ A). The resulting expression will represent a correct expansion 
of a K-dependent function in powers and logarithms of n to order o(/€^^). 



Concerning cancellations of the IR divergences in the first sum on the r.h.s. of [7.40| 



by the apparently UV divergent expressions in |7.41| , see the end of subsect. ^ 



Note also that |7.40| is a true infinite asymptotic expansion of (p^ + k^) ^ in the 



sense of distributions, and after the regularization is taken off, one only has powers and 
logarithms of k on the r.h.s. of |7.40|, and such an expansion is unique (cf. subsect. ^.2|). 



7.10 Several /t-dependent factors 

Let us consider an example with several factors, such that the singularities of their ex- 
pansions are located at different points. For example, consider the following integrand 
corresponding to a self-energy contribution: 

X 77 ^ ^. 7.42 



We regard the external momentum Q as a heavy parameter and fix it at a non-zero 
value, which corresponds to the simplified version of the expansion problem without 



contact terms (subsect. |3.5| ). The final recipe for expanding [7.42| in the sense of distri- 
butions will be to simply multiply the expansions for the two factors. But the reasoning 
behind it is more instructive than the result itself (cf. a similar reasoning in a more 
complicated case below in subsects. [13. 2| and |13.3|). 



The key observation here is that the expansion of the second factor is regular around 
the singular point of the first factor; more precisely, if ipi is such that it is equal to zero 
around the singular part of the second factor, then the ^-dependent test function 



satisfies the condition (ii) of subsect. |5.4| . Therefore, one can formally expand (pi and 
use |7.40| for (p^ + k^)~^, and we conclude that for p 7^ Q, the expansion of [r.42| is 



given by the Taylor expansion of the second factor times the full expansion of the first 
factor. Swapping the factors and considering test functions which are identically zero 
around the singular point of the first factor, we arrive at a similar conclusion. The final 
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observation is that since the singular points of the two factors are separated, any if{p) 
can be represented as fi{p) + ^2{p), whence the desired result can be easily deduced. 
(Another way to see it is to notice that products of ^-functions from expansions of the 
two products result in zeros.) 



To summarize, we have demonstrated with a simple example how the ideas of the 
extension principle discussed in the preceding section work in practice. We have also 
demonstrated a short-cut way to derive our final result, the expansion [7.40| , represented 
in the dimensionally regularized form. It can be represented in an explicitly finite form 
with the regularization taken off (for more on this see |^), but the form [7.40 
convenient for our purposes. 



IS more 



8 Applications to one-loop integrals 

The results of the previous sections can now be applied to one-loop integrals. We first 
consider the UV-convergent case, and then show how the case of UV-divergent MS- 
renormalized integrals can be reduced, essentially, to the convergent case. The most 
important non-technical element of our reasoning is a new definition of the MS-scheme 
in terms of subtractions of UV-asymptotic part directly from the integrand, prior to 
any momentum integrations (subsect. and pT[1 ). 



8.1 UV-convergent one- loop integral 

Consider the following UV-convergent one-loop integral: 

1 1 



One would expect (we discuss this in more detail below) that although the "test func- 
tion" (p2 -|- M^)"^ is not a true Schwartz test function, its decrease at p — oo is 
sufficiently rapid to allow one to substitute the expansion |7.40| for (p^ + k^)"^ into |0 



in order to get the asymptotic expansion for k ^ of the integral |8]^ as a whole. 
Doing this and integrating the 5-functions, one obtains: 



eq.8.1 = hm^^ j d^'p-x _ / rf^p _ x , (8.2) 



D^^u p' {p^ + M^y J p" {p^ + M^y 



1 1 

+ K 

The validity of |8.2| can be easily checked by explicit integrations. 



In order to make connection with what follows, let us briefly discuss why the rapid 
decrease of the "test function" (p^ + M^)~^ in ^]T| allows one to substitute |7.40| directly 
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into |8.1|. By definition, integrations over infinite regions involve a limiting procedure. 



Let us make it explicit in |8.1| . Take a smooth function such that 

= 1, for \p\ < 1, (8.4) 
0, for IpI > 2. (8.5) 



We can introduce the cut-off $(p/A) into 



1 [ ^p/A) 



eq.8.1 = lim / d^p — 1 ' . } . (8.6) 

^ A-ooi ^ p^ + K^\{p'^ + M^f] ^ ' 



Now it is fully correct to use |7.40| in order to expand the expression under the limit 
sign to order, say, o(k^). Then one has to check that (i) the limit A ^ oo exists for each 
term of the expansion of [j?' + /t^)~^ and [ii) the remainder which was o(/t^) prior to 
taking limit will be o(k^) after it. While ii) is obvious, to check iii) one can recall the 
representation |7.31| of the remainder term. The desired result easily follows therefrom; 
moreover, the proof here closely repeats the reasoning of the paragraph after |7.4 

The conclusion is that for absolutely convergent 1-loop integrals the asymptotic 
expansions are normally obtained by substituting the expansions in the sense of the 
distributions into the integrand, and then performing termwise integrations over the 
loop momentum. 



8.2 Separation of UV and IR divergences 

A remark is in order here. The expansion p.2| contains two terms that diverge as D ^ 4 : 
the second integral is logarithmically divergent for p — 0, and the third one, for p — > oo. 
However, by construction, the expansion as a whole is finite as D — > 4. Therefore, the 
divergences must cancel. But the peculiar feature of |8.2| is that one divergence occurs 



in a K-dependent integral, while the other one in an integral depending on M but not 
on K. This is not satisfactory because in the final respect in applications one has to 
deal with finite functions of k and M. 

As we have already noted, all our expansions can be represented in an explicitly 
convergent form which makes no use of the dimensional regularization, but this requires 
a special techniques (see pO[). However, it is still possible to rearrange |8]^ in terms of 
K- and M-dependent contributions which are separately finite, by using only notions 
known to the practitioners of perturbative QFT who are familiar with the techniques 
of UV renormalization within the MS-scheme. The price to be paid for this is that 
IR divergences will be canceled by UV-counterterms, but we have already explained in 
subsect. 17^ why such things can happen. 



Look at the third integral on the r.h.s. of p.2| . It is an UV divergent 1-loop Feynman 
integral which can be made finite by adding an UV-counterterm which can be evaluated 
e.g. in the MS-scheme. Such UV-counterterms are polynomials in masses and external 
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momenta, and in our case can be taken in the form: 



9 const n A 



i.7) 



where /i is introduced to preserve dimensionahty (cf. p.l5| below). Now, using 3]7, the 



expansion |8]^ can be identically rewritten as: 



eq.8.1 



lim-^ 

_D-4 



K 



+ 



M4 



+ 



const X /i 



D-4 



K^const X fi^~^ 



4) 



+ 



D-A 



+ 



(8.9) 



(8.10) 



p2 + 

Since the first and third square-bracketed terms as well as the expression as a whole are 
finite in the limit D A, so finite must be the second term too, which is not difficult 
to verify by explicit calculations. 

Note that by rescaling the integration variables, one can exhibit the dependence of 
various terms in |8.8| on k, and M : 

1 



eq.8.8 



1 



d^'p- 



1 



X 



d^'p- 



X 



(p2 + ly 

const X (/x/M)^-^ 



p" 



[d^p- 
J 



1 



(p2 + ly 

const X (yu//t 



D-A 



.D~4-\ 



0[K 



(8.11) 



f8.12) 



(8.13) 



p2 + 1 D-A 

We stress once more that the numerical connection of the UV counterterm in the 
third term and the counterterm in the second term, whose role is to cancel the IR 
divergence — this connection can be traced back to the fact that the dimensional reg- 
ularization nullifies certain integrals containing both IR and UV divergences simul- 
taneously. One could pursue the topic further and analyze such cancellations more 
explicitly, but as the relevant formulae are quite complicated while the results per- 
fectly useless, we stop here. 

To conclude, we have demonstrated a trick for recasting expansions that emerge 
from our methods into an explicitly convergent form. The trick employs specific prop- 
erties of UV renormalization of Feynman integrals, and can be generalized to the 
multiloop case where it takes the form of an inversion of the i?-operation [p!^- [|l^ . 



8.3 UV divergent one-loop integrals 

Let us now turn to UV divergent one-loop integrals, e.g. 

1 1 



X 



(8.14) 
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Its MS-renormalized version which we wish to expand in k, is: 

1 1 z 



hm 



p2 + ^2 p2^M^ D-A 



D-4 



(8.15) 



where 2; is a numerical constant and fi is 't Hooft's unit of mass which plays the role 
of the subtraction parameter in the MS-scheme. 

Our purpose is to show that the expansion [7.27] (and the more general expansion 
|7.40|) can be directly substituted into the integral in |8.15| in order to get the asymp- 
totic expansion for ^.15| . (Correctness of this recipe can also be checked by explicit 
calculations.) The motivation comes from the fact that the divergent contribution to 
the integral in ^.15| is generated by the asymptotics of the integrand at p 00, while 
replacing [p^ + k"^)^^ by its expansion leaves that asymptotics intact. These heuristic 
arguments can be formalized in the following way. 

The crucial step is to redefine the MS-subtraction in |8.15| in a manner independent of 
dimensional regularization (because the expression B.15| does not allow one to apply the 
reasoning like that in subsect. |8.1| ). To this end recall (subsect. [7.7|) that a subtraction 
of large-p asymptotics from the integrand is equivalent to UV renormalization. One 
has: 

11 1 , , 
X ~ f8 16) 

^ ^ p 00 ^ 

{k,M 0) 

so that subtracting from the integrand in ^.15| would produce an UV convergent 
integral. However, p~^ is ill-defined at p = 0. Let us modify it at p = so as to ensure 
its integrability: 



p4 



R 



■fx 



.D-A 



Sip) 



(8.17) 



p4 D-A 

(cf. [7.10 ). z' in p. 17] is a numeric constant, and the mass-like parameter /i is introduced 
to preserve dimensionality. 

Now, the expression 



is finite at D 



d^p 

A both at p — > 
d'^p 



P^ + 

and p - 
1 

X 



1 


" 1 " 




p2 + M2 




"1 



00. Owing to |7.29| , it can be represented as: 
1 z' 



P^ + 



p^ + M"^ D-A 



D-A 



(8.19) 



Recall that the constant z in p.l5| was chosen solely from the requirement of finiteness 
of |8.15| . Hence, comparing j.l5 and 8.19 we see that z = z\ and eq. jj.l5 = eq. 8.18| . 



Finally, recall that the distribution |8.17| can be represented as |7.8| with a suitably 
chosen c, and no regularization is used in [7.8|. Therefore: 



eq.8.15 



d^p 



1 



p4 



(8.20) 



8 APPLICATIONS TO ONE-LOOP INTEGRALS 



45 



So, we have represented the MS-renormahzed integral ^.151 in the form of an integral 
over p of the expression which differs from the original integrand ^.141 by a contribution 
independent of k. 



Further reasoning presents no problems, as the representation |8.2CI| allows one to 
use the arguments of subsect. with only inessential modifications in order to obtain 
the expansion in n. We leave it to an interested reader to check that the final result is 
such as if we had substituted the expansion flM directly into |8J[5| . 



The method for recasting the expansion into an explicitly convergent form described 
in subsect. is also applicable here. 



8.4 Generalized MS-schemes 

Now consider the following definition of a class of subtractions schemes. (It is shown 



elsewhere that this definition works in the case of multiloop Feynman integrals ^1 
|1T9[| .) Let F{p) be the integrand of a 1-loop Feynman diagram. Evaluate its asymptotic 
expansion for p ^ oo and discard the terms that are integrable at large p. Denote the 
result as F^^{p). The terms in F^^{p) that generate logarithmic divergences at p — cxd, 
also have a logarithmic singularity at p = 0, while the rest are integrable at that 
point. Redefine the logarithmically divergent terms at p = similarly to [7.8| , [7.101 
to make them integrable at all finite p. Denote the result as [F'^^{p)]^. (Note that 
the redefinition should only affect powers of p and commute with other dimensional 
parameters; in particular, the constant c in |7]^ should be independent of the masses 



and external momenta of the diagram.) The renormalized integral is now defined as 

/ d'p {Fip) - \F-^{p)f) . (8.2i; 



It is easy to see (cf. the reasoning in subsect. ^.3[ ) that the class of subtraction schemes 
thus defined includes the MS-scheme as a special case. Note also that the UV- 
counterterms in these schemes are always polynomials in masses and external momenta 
of the diagram. 

To take an example, consider the integral that has already been encountered above 
in subsect. fTT] : 



then 



and 



p2 _|_ ^2 ' p2 pi ' 



p4 



(8.24) 



The renormalized integral is 
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Introducing the dimensional regularization into |8.25| , using p.l7| and |7.29| , we arrive at 
an expression, 

where z' = Cq + {D — 4)ci + . . . and q are arbitrary finite constants. If q = 0, i > 1, 
one recovers the standard MS-scheme. 



8.5 Summary 

The representation of the MS-renormalization as subtractions from the bare integrand 
of those terms from its asymptotics at large loop momentum which are responsible 
for UV divergences, allows one to reduce the problem of expansion of 1-loop integrals 
to the UV convergent case. And the reasoning of the above sections shows that the 
correct expansion of the renormalized integrals can be obtained by direct substitution 
of the corresponding expansion in the sense of distributions in place of the K-dependent 
terms, and rearranging the resulting expression along the lines of subsect. p.2| , in order 
to obtain the expansion in an explicitly convergent form. 

In order to generalize this recipe to the multiloop case we need, first, to learn to 
expand more complicated products of singular functions (which is done in the present 
paper), second, to generalize the above treatment of UV-divergent integrals to the 



multiloop case [|2T|, [|T^], and third, to generalize the rearrangement trick of subsect. |8]2 



THE MASTER EXPANSION 



9 Notations for products of singular functions 

Admittedly, multiloop Feynman diagrams are cumbersome objects. Nevertheless, the 
analytical aspects of the problem, say, of UV renormalization or the Euclidean expan- 
sion problem are by no means as complicated as one can imagine regarding them from 
the point of view of the traditional techniques prevailing in rigorous studies of per- 
turbative QFT like e.g. the a-representation which completely destroys the essential 
multiplicative structures of Feynman integrands. However, to reveal the underlying 
simplicity and help one to ignore the plethora of irrelevant detail, a careful choice of 
notations is crucial. Below we describe such a system of notations (discussed in a more 
formalized manner in [^) which does not make use of parametric representations and 
applies equally well to both scalar and non-scalar theories. 
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9.1 Graphs and products 

Let G be a Feynman graph. Its loop momenta pi,p2, ■ ■ ■ ,pi are always Euclidean vectors 
of 4 dimensions. (For studying toy models and examples it is sometimes convenient to 
consider the 2-dimensional case.) The loop momenta are combined into a single vector 
of 4 X Z dimensions: 

P = iPl,P2- ■ -Pi) ■ (9.1) 
The variable p may also incorporate the heavy external momenta — see remark {in), 



subsect. In such a case 

P= {pi,P2---PhQi---) , (9.2) 

where Q are defined in subsect. ^]4|. 

The integrand of the graph G in momentum representation is denoted as G{p). 
Dependence on other parameters can be shown with additional arguments after p. 
G{p) is a product of factors corresponding to both lines and vertices. The factors are 
numerated by the label g. The set of values of the label is denoted also as G. The use 
of the same symbol for three different but related entities should not be misleading. 
Similarly, it should not be misleading if we denote as g{p) that factor in the product 
G{p), which corresponds to the label g. So, 

G{p) = n g{p). (9.3) 

geG 

We will also have to study subproducts of G. Let 7 C G. Then 

lip)^Yl9ip) (9-4) 

and also 

G\7(p) = n 9ip) = G{p)h{p). (9.5) 
where G\7 is the standard notation for the difference of sets. 



9.2 Conventions for non-scalar factors 

We allow each function g{p) to carry sub- (or super-) scripts: Lorentzian, SU(?t,) etc. 
Then the product |9^ may contain implicit contractions. If, for example, the product 



carries a pair of subscripts, which are contracted in p.3| , then this contraction is 
implied in |9.4| as well. 



Now, let F' and F" be non-intersecting subproducts from G. In G there may exist 
such a pair of contracted subscripts that one of them belongs to F', while the other 
one to F". Then the product 

V\p)V"{p) (9.6) 
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implies contraction of such a pair. Furthermore, we will build new expressions from G 
by replacing some subproducts 7 by certain expressions, say, Z^, so that will carry 
exactly the same subscripts as 7. Then the products like 

r'Zr", Zr'Zr" etc. (9.7) 



imply the same contractions as pl3 



The implicit presence of such contractions will have no complicating effect on our 
reasoning as compared with the scalar case. 



9.3 Momentum dependence 

Let us discuss the dependence of factors on p. One should distinguish factors of two 
types: those corresponding to lines and to vertices of Feynman graphs. 

Let g correspond to a line. The momentum flowing through this line is a linear 
combination of various momentum variables of our problem: 

( E Cg,iP^ + Qg] +kg, (9.8) 
\i=l...l / 

where Qg and kg are linear combinations of heavy and light external momenta from the 
sets Q and k. Then g{p) is a propagator depending on and a mass. The mass may 
be light, i.e. equal to zero or proportional to k, or heavy, i.e. non-zero and independent 
of K. Since kg in |9.8| is proportional to k, and will in fact play the role of an expansion 



parameter, it is convenient to introduce a notation for the non- vanishing part of pl8: 
Lg{p) = the parenthesised expression in 9.8. (9.9) 

Now one can write: 

gip)=FgiLgip)+kg,K). (9.10) 

In fact our arguments are true for a wide class of functions Fg{L,K) including any 
standard perturbation theory propagators raised to integer powers and multiplied by 
polynomials of L. 

9.4 Assumptions on the properties of factors 

Let us describe those properties of Fg{L,K,) that are essential for our theory. (For a 
formalized description see [^.) One should consider two cases: 



(i) The function Fg{L, k) at L = k = is regular (which corresponds to propagators 
with heavy masses). Then for k 0, Fg is expanded into a series of functions which 
are smooth in L, and the expansion is uniform in L in any bounded region and allows 
arbitrary termwise differentiations in L. 
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(ii) The function Fg{L, k) at L = k = is singular (which corresponds to propaga- 
tors with a hght mass). Here we assume that Fg{L, k) has the following properties: 

(a) scaling: 

Fg{XL, \k) = X'^'FgiL, k); (9.11) 

(b) asymptotic expansion at each L ^ : 

Fg{L, k)=Y. K^Fg,^{L), L^O; (9.12) 

n 

(c) all Fg^n{L) are smooth at L 7^ and have simple scaling properties: 

Fg,„(AL) =A'^''-"Fg,„(L); (9.13) 

(d) regularity properties: the expansion in (b) allows arbitrary termwise differenti- 
ations in L, and the remainder of the expansion to order 0{n^) has an upper bound 
of the form of the modulo of the first discarded term: k^^^\L\'^o~^~^ (cf. property (c); 
also note that property (d) ensures that the ^-dependent test functions that arise at 
intermediate steps of our reasoning will possess the properties described in subsect. p.4p . 

In case i/ii we will say that the factor is of regular / singular type, respectively. 

It should be stressed that the described properties are sufficient for our purposes 
and comprise the problems within the standard perturbation theory. But our methods 
will work in a more general context, e.g. if the singular functions have logarithmic or 
power behaviour with a non- integer exponent. 



9.5 Factors corresponding to vertices 

Now, let g correspond to a vertex of the graph. One can assume that g{p) is a uniform 
polynomial of momenta entering the vertex. There are three ways to include such fac- 
tors into a common framework with lines. First, all vertex polynomials can be merged 
with test functions. Second, each vertex polynomial can be split so that each part may 
be included into one of the propagators attached to the vertex (such modifications of 
propagators are allowed by our definitions — see subsect. |9^) . The third way — which 
we prefer — consists in generalizing the notations introduced for lines: it is sufficient to 



allow Lg to consist of several independent momentum components like 9.8, i.e.: 



^9=\Y^ Cg,v,iPi + Qg,v\ (9.14) 
\'i=l...l / t)=l 

For convenience, we assume that, by definition, the vertex functions belong to the 
singular type (see subsect. above). 

Equivalence of expansions obtained in all three cases follows from the uniqueness 
of the final result (subsect. |5^) . More explicitly, the equivalence can be checked using 
the property of the As-operation to commute with multiplications of the expression to 
be expanded by polynomials of p (see below subsect. p.2.5|) . 
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9.6 Summary 



Let us summarize the entire scheme. The product p.3| is defined on the space of the 
multicomponent variable |9.14| , and consists of factors of the form |9.10| . The factors 
need not be scalar functions — the conventions of subsect. are operative. The first 



argument in |9.10| is a set of momentum variables of the form |9.14| . The functions Fg 
satisfy the conditions described in subsect. |9]^. Finally, the factors may be of two 
types — singular or regular (see the end of subsect. [9^). 



10 Formal expansions and IR-subgraphs 

An important element of our technique is classification of singularities of the formal 
(Taylor) expansions in k of the products |9.3| . The key notion here is that of IR- 
subgraph, introduced in subsect. |10.4 Note that there are similarities between our 
Euclidean space classification of IR-singular points and the Minkowski space reasoning 



of Libby and Sterman |66 



10.1 Some notations 

Denote the operation of formal expansion in k, as T^. The expansion Tf^g is the same as 



the one described in property (b), subsect. up to a replacement of the momentum 
variable. Then 

T«^7(p,«:)= Y[T,gip,K), (10.1) 

and similarly for each subset 'j G G. The r.h.s. of |10.1| should be understood as a 
simple infinite series in integer powers of k, obtained by formal multiplication of the 
series with reordering of terms. 

It is convenient to denote the operation of partial expansion up to terms O(k^) as 
T^K,N- We will also omit the subscript indicating the expansion parameter as —>■ T, 
AT — >■ Ttv, if misunderstanding is excluded. 



10.2 Singularities of the formal expansion 

Let us describe the structure of singularities of the formal expansion |10.1| . Each factor 
of the singular type generates singularities at those points where Lg{p) = 0. We call 
the set of such points singular plane of g and denote it as 

7,g = {p\Lgip)=0}. (10.2) 

Now consider an arbitrary subset 7 C G such that all its elements are of singular type. 
Define the singular plane of the subset 7 : 

vr, = n rrg. (10.3) 

367 
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On this plane, all factors (767 are singular simultaneously. 

We will consider only such subsets 7, to which there correspond non-empty singular 
planes. 

An important point is that the same singular plane may correspond to different 
subsets. For example, in Fig. 3 the pairs of propagators with momenta pi and p2', Pi 
and pi + P2 generate the same singular plane described by the equations Pi = P2 = 0. 
Further, there is one largest among the subsets with the same singular plane. This 
follows from the fact that if vr^ = n^i then In Fig. 3 such a subproduct 

consists of all three propagators. Subsets, to which new factors cannot be added 
without reducing their singular plane, will be described as complete and will be called 
(complete) IR-subgraphs or, if confusion is excluded, simply subgraphs. 

10.3 Diagrammatic interpretation of IR-subgraphs 

Take a Feynman graph and evaluate momenta for each line. To check if a given set of 
singular lines 7 is complete, perform the following test. Set to zero all the momenta 
flowing through the lines of 7, and all light external momenta. Using momentum 
conservation at vertices, reevaluate the momenta for the rest of the lines. If there is a 
line of singular type not belonging to 7, whose momentum will vanish after reevaluation, 
then the set 7 is not complete. If there are no such lines, then 7 is an IR-subgraph. 
For example, in Fig. 4a the pair of vertical lines forms an IR-subgraph. 

Note that the same subset 7 in the same diagram G may be or be not an IR- 
subgraph, depending on whether the full Euclidean expansion problem is considered 
or its simplified version without contact terms (see subsect. |3.5| ). For example, for the 
graph in Fig. 5 (masses in both lines are light) the entire graph will be its own IR- 
subgraph if the graph is considered as a distribution over Q, so that Q is included into 
integration momenta; the corresponding singular plane is described by the equations 
p = Q = 0. In the simplified version one considers integration only over the loop 
momentum p, while Q is fixed at a non-exceptional value (this means here simply that 
(5 7^ 0); in this case the pair of lines cannot be an IR-subgraph because its singular 
plane is empty. 

The property of being IR-subgraph also depends on the additional restrictions spec- 
ified for the heavy momenta in the particular problem. For example, in Fig. 4a which 
differs from Fig. 4b by a linear restriction Qi = Q2, the pair of vertical lines forms an 
IR-subgraph, while this is not so for Fig. 4b. 

When should a vertex be included into an IR-subgraph? Recall that any vertex is a 
singular-type factor by definition (cf. subsect. and its singularities are localized at 
the points where all the momenta from the set Lg are equal to zero. Therefore, in the 
most general case the recipe is as follows: if all the lines incident to this vertex belong 
to the subgraph, then the vertex is to be included, too. However, in the simplified 
version of the expansion problem, an IR-subgraph may not contain a vertex with a 
non-zero total heavy external momentum entering into it. 
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Starting from the basic analytical definition, one can always enumerate IR-sub- 
graplis in any specific case. However, we do not attempt translating the definition of 
the IR-subgraph into the graph-theoretic language: the result would be cumbersome 
and practically useless. For the purposes of doing the combinatorics of factorization 
1T8| , we will only need the properties of complements of IR-subgraphs that can be 



m 



easily derived from the above definition (cf. subsects. f).2[ |6.1| and |6.4| in |jT 



10.4 Why IR-subgraphs? 

The reason for considering only IR-subgraphs among all subproducts of G(j)) is as 
follows. If 7 is a subgraph, then singularities of the formal expansion 

T,G{p) = [T^^ip)] X [T.G\7(p)] (10.4) 

near almost any point at the plane vr^ (except for the intersections of with singular 
planes of factors not belonging to 7 ), are determined solely by 7. 

Vice versa, let po be any point from the space of p. Select from the graph the 
smallest subproduct 7 such that TkG\7(p) have no singularities at p = pq. Then 7 will 
be an IR-subgraph. 



10.5 Hierarchy of IR-subgraphs 

The set of all IR-subgraphs of a given graph has a natural ordering. Namely, 7 < F 
if and only if 7 belongs to F. Note that if 7 < F, then vr^ D vrp. Correspondingly, one 
can define maximal and minimal subgraphs. Note that there may be several maximal 
subgraphs; their singular planes do not intersect. 

It is convenient to introduce the notion of immediate precedence for subgraphs. We 
say that a subgraph 7 immediately precedes the subgraph F, and write 7<F, if 7 < F 
and there are no subgraphs between 7 and F, i.e. no subgraphs F' such that 7 < F' < F. 
We describe the subgraph which immediately precedes a maximal one as submaximal. 



10.6 Proper variables of subgraphs 

Let F be a subgraph of the graph G. In general, T{p, n) is independent of some com- 
ponents of p. Indeed, it is easy to check that T{p, k) is always invariant with respect 
to shifts in p parallel to vrr. Split p into two parts: transverse and longitudinal with 
respect to vrr : p = {j)v,Pr)i so that the singular plane vrr is described by the equation 
Pr = 0. Then F(p, k) depends only on p^. The variables pr will be referred to as the 
set of proper variables of F. 



Although the choice of proper variables is not unique, transition to another set 
of proper variables is equivalent to a change of coordinates in the space of pr- Here 
invariant coordinateless formulations are possible and useful — for more on this see . 
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The proper variables of an IR-subgraph can be easily determined using its graphical 
representation. Thus, for the subgraph in Fig. 3 the proper variables are the set of 
momenta pi and p2 (or pi and P2 = Pi + P2 etc.). 

Consider the case when the IR-subgraph has loops (see e.g. Fig. 6). The loop 
momenta belong to the proper variables of the subgraph. Denote the set of all loop 
momenta as p™*, and the rest of the proper variables as p^^^*. So, 

P,= {p'^\pT)- (10-5) 



To conclude this section, note that any IR-subgraph can be considered as a graph in 
its own right. Then the loop momenta of the graph should comprise all proper variables 
of the IR-subgraph. Graphically, this corresponds to merging all the "external" vertices 
of the IR-subgraph into one. (In Fig. 6, the external vertices are a and b.) The resulting 
graph possesses an additional property: it is its own (and the only maximal) IR- 
subgraph. In that case the singular plane of the entire graph is reduced to a point which 
can be assumed to be the zero point. It is clear that for each subgraph one can consider 



the problem of its asymptotic expansion. Moreover, we will see in subsect. p.l.3| that 
construction of the As-operation for the graph is logically preceded by its construction 
for the subgraphs. It is worth noting that the As-operation for a subgraph is formulated 
entirely in terms of its proper variables. 



10.7 Quantitative characteristics 

Let us introduce some quantitative characteristics for description of singularities of 
IR-subgraphs. For each IR-subgraph F define: 

dr = J2d9, (10.6) 

which describes simultaneous scaling in k, and the proper variables: 

r(Ap,AK) = A*r(pr,/«) (10.7) 

(cf. property (a), subsect. |9.4D; 



dimpr (10.8) 
is the dimension of the space of proper variables of the subgraph; 

ur = —dr — dimpr (10.9) 

is the singularity index. Its meaning is as follows. If for all F < G one has uJr + 

< (strict inequality!) then the formal expansion of G to order does not 
contain non-integrable singularities and is a correct asymptotic expansion in the sense 
of distributions. 
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11 As-operation for products of singular functions 

In this section, following the recipe of the extension principle (sect. 5), we establish 
formulae for the expansion of products of singular functions in the sense of distribu- 
tions in the form of the so-called As-operation whose structure is similar to that of 
the Bogoliubov i?-operation. The derivation presented bears a resemblance to, and 
was partially inspired by the analysis of the i?-operation in |2^. Explicit expressions 
for counterterms of the As-operation (which, unlike the i?-operation, is determined 
uniquely) will be obtained in the next section. 



11.1 General remarks and notations 

We are going to derive the recipe of the As-operation, the instrument for evaluating 
expansions in the form of infinite asymptotic series for any product (graph) G{p) of 
the described type. We will use the notation: 

G{p,^) - AsG(p,«:) =^«:"G„(p,«:). (11.1) 

Here each Gn should be a distribution in p, and the dependence on n should be soft 
(see the definition in subsect. p.lj) . The expansion |11.1| should be valid in the sense of 
distributions (see subsect. p.l| ). This means that for each test function ip the following 
estimate should hold: 



J dp^ip) [G{p, k) - As^Gip, k)] = o(/s:^). (11.2) 



The subscript on the As-operation means that in |11.2| all the terms of order o(k^) are 
discarded, i.e. 



AsatG = J2 ^"Gn- (11-3) 

n<N 

Existence of such expansions is either obvious or non-trivial, depending on the point 
of view. It is rather obvious if one recalls the results of [|l| wherefrom one concludes 
that expansion of G integrated with a test function runs in powers and logarithms 
of the expansion parameter and takes into account that such an expansion should be 
linear in test functions. On the other hand, from the point of view of the technique 
of PI existence of such expansions is not at all trivial because the construction of H] 



corresponds to an expansion similar to |11.3| for each but with Gn depending on A^, 
too. 



11.2 Uniqueness of the expansion 



From subsect. ^]2| it follows that if an expansion of the form |11.1| exists, then it is 
unique (provided one always expands in powers and logarithms of the small parameter). 



Uniqueness of the As-operation will allow us in subsect. |11.3| to determine the necessary 
conditions which it must satisfy. This will give us a sufficient number of hints to 
construct it explicitly. 
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11.3 Localization property of the 74s-operation 



Let us study the local structure of the As-operation (cf. remark (i) in subsect. O). 
First of all note that the formal expansion 



:il.4) 



allows integrations with test functions ip whose support does not intersect singular 
planes of the formal expansion: 



suppy^nfU^ 



;ii.5) 



This follows from the properties of the functions of which G is built (cf. subsect. |9.4|) . 
The r.h.s. of 11.4 defines a functional which is an asymptotic expansion of the functional 
on the l.h.s., on the subspace of test functions satisfying |11.5|. From uniqueness it 



follows that on test functions satisfying |11.5| the As-operation should coincide with the 
operation of the formal expansion T: 



AsG{p,k) = TG{p,k) 



;ii.6) 



Further, consider any small region O, and let F be a subgraph such that TG\F(j9) 
has no singularities within O (cf. subsect. [4.1|) . Take an arbitrary ip such that supp C 
O. Then 

I dpG{p,KMp) = J T{p,K)ip"{p,K), (11.7) 

where 

if"ip,K) = G\Tip,KMp) (11.8) 

is a test function with supp ip C supp ip", but depending on k. 

Recall (subsect. |10.6|) that F is independent of that component of p which is lon- 
gitudinal with respect to vrr. Therefore we may perform integration of ip" over the 



longitudinal components on the r.h.s. of 11.7. So, we arrive at the problem of expan 



sion for the subgraph in terms of its proper variables, and with test functions depending 
on the expansion parameter. The assumptions made in subsect. |9.4| ensure that the 
expansion procedure described in subsect. |5]J will be applicable here. 



So, assume that the problem of expanding F(p, k) is solved, i.e. that the action of 
the As-operation on F is already known. Then the expansion of the r.h.s. of |1 1 . 7| can 
be represented as 

J dp [AsF(p, k)] X [Tif"{p, k)] = j dp [AsF(p, k)] x [TG\T{p, k)] ^{p). (11.9) 



On the other hand, the expansion of the r.h.s. of |11.7| should be given by 



J dp As G{p,K)ip{p). (11.10) 
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Then again from uniqueness of the As-operation we conclude that, given the relation 
described above between the region O and the subgraph T, on all test functions localized 
in O the following localization property must be valid: 

AsG{p) = [TG\r{p)] X [Asr(p)], peO. (11.11) 

This is the most important structural property of the As-operation. It exhibits the 
recursive structure of the expansion problem considered on the entire collection of 
Feynman diagrams. 



11.4 Structure of the 74s-operation 

Let us represent the As-operation in terms of counterterms localized at singular points 
of the formal expansion, in analogy with the well-known expression of the Bogoliubov 
i?-operation in terms of quasilocal counterterms p4|. AsG may differ from TG only 
by corrections localized on singular planes of the formal expansion. Consequently, one 
can take the following ansatz for the full expansion: 

AsG = TG + ^(E7)irG,7, (11.12) 

7 

where summation runs over all IR-subgraphs, while the distribution (E7) is localized 
on TT^. Let us determine K^-y from |11.11| , expanding the As-operation in accordance 
with ITTT^ : 

TG+Y: (E7) Kg,, = (TG\r) (TF + ^ (E7) i^r.^)- (11-13) 

7<G 7<r 

We get the conditions: 

XG,7 = irr,7(TG\F). (11.14) 
Without loss of generality we may assume that for each F : 

Kr,r = 1. (11.15) 

(Indeed, within the problem for the subgraph F considered in terms of its proper 
variables, the distribution (EF) is localized at the origin, so that K^^r must be a 
constant independent of momentum variables. And such a constant can always be 
included into (EF).) 

Setting 7 = F in |11.14| and using |11.15| , we get: 

KG,r = TG\T (11.16) 

for each subgraph F < G. We see that the dependence on F disappears from the r.h.s. 
of |11.14| , as expected. Finally: 

AsG = TG+ ^ (E7)(TG'\7). (11.17) 

7<G 



The expressions (E7) will be referred to as counterterms for subgraphs 7. 
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For |11.1| to be true, it is necessary that the counterterm had the form of an infinite 
series 

(E7) = 5: (E7). , (11.18) 

n 

where (£7)^ depend softly on k. 



11.5 Summary 

So, if the As-operation exists, it is unique and, in accordance with the exten- 

sion principle, must have the form |11.17Hri.l8| . However, the analysis of subsect. |11.5 



indicates that the As-operation can be naturally defined by [11 .171 using the induction 



based on the natural order among IR-subgraphs (see also below sect. |T2|). Then if the 
As-operation has been constructed for all 7 < F, it only remains to determine the 
counterterms (EF) possessing all the required properties. Note that the reasoning of 
sect. is in fact applicable to any minimal IR-subgraphs, which provides a correct 
starting point for the induction. The general formulae for the counterterms derived in 
the next section will include the case of minimal subgraphs as a simple special case. 



12 Counterterms of the As-operation 
12.1 Structure of the As'-operation 

Suppose that the existence of the As-operation has been established for all 7 < F, and 
explicit expressions for (E7) of the form of |11.18| have been found. From the reasoning 



of the preceding section it follows that the expansion for F valid on test functions which 
are equal to zero around p = has the form: 

T{pr,K) - As'T{pr,K), pr ^ 0, (12.1) 

where As' is the As-operation without the last counterterm: 

As'F = TF+5:(E7)(TF\7), (12.2) 

7<r 

so that 

AsF = As'F + (EF). (12.3) 
Indeed, let H run over all submaximal subgraphs of F, i.e. H <T (see subsect. |10.5|) . 



Then pairwise intersections of the singular planes tth consist only of the point pr = 0. 
Therefore, if ip{pr) = around pr = then ip can be represented as 



H<\r 



(12.4) 
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where <fH{pr) = in a neighbourhood of tch' for any H' ^ H. Then 

<Vip> =Y.<T^H>=J2<H,T\H^H> 
H<ir H<ir 

J2 < As H,TT\HipH> 



;i2.5) 



;i2.6) 



(cf. the reasoning in subsect. |11.3| , especially |11.9|) . On the other hand, the operation 
As' as defined in p.2.2| inherits the localization property analogous to 



< AsT,ipH >=< AsH, TT\HipH >, 



so that the r.h.s. of 12.5 is equal to 



< AsT, (fH >=< AsT, (p >, 

H<ir 



[12.7) 



:i2.8) 



whence follows 12.1 



12.2 Approximation properties of the As'-operation 

Our inductive assumptions imply that all counterterms (E7) for 7 < F have the form 



11.18| . Then As'F can be represented as an expansion in powers of k : 

As'r(pr,/«)=E'^X(Pr,«:), 



:i2.9) 



where T'^{pr,K.) can softly depend on k. Besides, the functions possess a number 
of natural properties following from our assumptions. First, they inherit the scaling 



property |10.7 



F;(Apr,A«:) = A'^^-"F;(pr,«: 



;i2.io) 



Second, from |12.10| and the soft dependence on k, one immediately gets the scaling 
property in the momentum argument: 



F^(Apr, = A'*" " (r^(pr, f^) + soft corrections) 



12.111 



To put it differently, the dependence of the expression A" »F^(Apr, 1^) on A is soft. 

Note the important parallel between what we have here and what we had in sub- 
sect. |7.1| : the expansion |12.9| , similarly to [7.2| , contains terms progressively more singu- 
lar at pr ~ as n —> 00, and both expansions are valid in the sense of the distribution 



theory for pr 7^ 0. We can push the analogy even further. Consider the expansion p.2.£ 
truncated at the terms of order O(k^) : 



F = As'^F + o{k 



(12.12) 



(The use of the subscript here is analogous to |11.3| .) From |12.11| it follows that the 
r.h.s. of |12.12| becomes non-integrable at zero for u-p + N > (cur is defined in |10.9|) . 
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However if the test function has a zero of order cur + iV + 1 at pr = 0, then the r.h.s. 
of |12.12| is well defined and, moreover, on such test functions the estimate o(ft^) of the 
expansion 



12.12 is valid. 



The last property is fully natural. Its proof is essentially not difficult, is entirely 
based on power counting and has the same structure as that in subsect. \I72^ One 
considers a neighborhood of the origin of a radius 0{k) wherein the uniformity prop- 
erty |12.10| is used for explicitly extracting the factor 0{k^~^^) from the integral of the 
remainder term of the expansion over this neighbourhood. Then one uses the prop- 
erties of the As-operation on subgraphs (the inductive assumption) to estimate the 
integral over the rest of the space. This last step is somewhat cumbersome though 
fully straightforward. A detailed discussion can be found in [^. (See also sect. |13 
below.) 



12.3 Expressions for counterterms 



The above results allow one to apply here the same reasoning as in the example in 
sect. 1^. From the extension principle it follows that the expansion |12.12| can be extended 
to all pr by addition of counterterms localized at pr = : 

r = As^^ + (E^.^) + o(/^^). (12.13) 

(E^vF) may only contain derivatives of the ^-function of order not higher than uj^ + N : 

(E^F) (pr) = E crM6r,aiPr) (12-14) 

(here 5r,o(pr) is a full basis of 5-functions and their derivatives localized at pr = 0; the 
order of derivatives is denoted as |q;|). To find explicit expressions of the coefficients 
cr,a (recall that in general they depend on A^, but it will turn out possible to choose 
them independent of A^), we use the trick of subsect. fTS . 

Consider the expression 



J dprA{pr)B{pr) = o(«:^). 



where 



A 
B -- 



F - As'^F, 



12.15) 

;i2.i6) 
;i2.i7) 



with C(Jr as defined in |1U.9| . The fact that the l.h.s. of |12.15| is indeed o{k ) is proved 
exactly as in subsect. |7.9| ; one only has to note that owing to the scaling properties 
TOTi , |l2l0| one has: 



(/ - As'^) F(Apr, = (/ - As^) F(pr, 



A— +00 



((/./A) 



N 



Pr^O. 



;i2.i8) 

12.19) 
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Therefore, the integral converges at pr oo by power counting. (Note, however, that 
this is not an ordinary absolutely convergent integral because the term As^F is a 



distribution. A more accurate definition of the integral in 12.15 involves introducing a 
cut-off $(pr / A) with any $ defined as in ^.4| and taking the limit A oo. The resulting 
value is independent of the exact form of Also note that the integral |12.15| can be 
transformed into an absolutely convergent integral of ordinary functions by invoking 
the explicit expressions for the operation As' and integrating out the 5-functions.) 



Introduce the dimensional regularization into the l.h.s. of P-2.15| , multiply A and B, 
and perform termwise integrations using the fact that the dimensionally regularized 
integration nullifies integrands with simple scaling properties — cf. 7.29| . Finally, we 
arrive at the following analogue of |7.41| : 



12.20) 



where Vr,a is the basis of polynomials of the variable pr, which is dual to the basis of 
5-functions 5r,a(pr) — cf. [7.37| . The counterterms (EF)„ from pTTTSl now are as follows: 



(Er)„ 



J2 cr,, 

\a\=LJi'+n 



r,a- 



12.21) 



Using the scaling properties, we get: 



crM = K(l"l-i'+-)cr,„(l), 



ri2.22) 



where the dots denote a correction proportional to D — 4, the deviation of the complex 
parameter of dimensionality from the canonical integer value. The integer part of the 
exponent in |12.22| shows which order in k, this term belongs to. One can also see that 



the dependence on k in the expansion AsF is as discussed in subsect. 5.3 



The formulae |1 1 . 1 7| , |11.18| , |12.21| and |12.20| present the full solution of the Master 
problem in the Euclidean case. 



12.4 Remarks 

(i) It should be stressed that although the above formulae are represented in a form 
which is specific to dimensional regularization, there are no serious obstacles (except for 
necessity to introduce a large number of new notations) for writing analogous formulae 
in other regularizations (e.g. using cut-offs). But this falls outside the scope of the 
present paper (see, however, [pO[] ). 

(a) Concerning the mathematical aspects of our derivation, the full details of the 
proof that are regularization independent are presented in pO|. As to the use of the di- 



mensional regularization in the preceding section, it can be justified along the following 
lines. First, the definition of dimensionally regularized integrals of arbitrary smooth 



test functions over momentum variables was presented in BB]. Second, as noted after 



12.18| , the integral in |12.15| can be represented as an absolutely convergent integral 
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of ordinary smooth functions. Therefore, one only has to extend the results of |35 



to a wider class of functions, which seems to present no difficulties. This issue will, 
hopefully, be addressed in a future publication |3^. Another approach would be to 



check the final formulae for expansions of Feynman diagrams using the a-parametric 
integral representation. But the attempts to do so show that a full proof would be 
a straightforward, cumbersome and unilluminating excercise. 



12.5 Commutativity with multiplication by polynomials 

In conclusion of the discussion of the Master problem, a remark is in order. We have 
proved existence of the As-operation and its uniqueness (subsect. ^.2|) . It follows that 
it commutes with multiplication of F by an arbitrary polynomial of momenta — indeed, 
the polynomial can either be included into the distribution to be expanded (F) or, 
equivalently, it can be considered to be a part of the test function; the final result cannot 
depend on which option is chosen. This means that the operation E which generates 
the counterterm for a subgraph, commutes with multiplication by polynomials, too. 
This can also be checked with the help of the explicit representations obtained above. 

The described property of the As-operation is useful in situations when one needs 
to transform expansions of Feynman diagrams to a convenient form. 



13 Example of an ^5-expansion 



13.1 The expression to be expanded and its singularities 

Let the integration variable p consist of two 4-dimensional components: 

P={Pl,P2)- (13.1) 

The product to be expanded is (cf. Fig. 7a): 

Gip, k) = gi{p, k) g2{p, n) g3{p, k), (13.2) 

where ^ 

9i{p,f^) = —, — 2^ (13-3) 

Pl + 
P2 + K. 

93{P, 1^) = 7 ^2^ 2 - (13.5) 

[Pl - P2Y + 

The expansion of an isolated propagator in the sense of the distribution theory has 
been obtained and studied in sect. ^. 
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Each propagator in |13.2| generates, upon formal Taylor expansion in k, singularities 

at the points where the corresponding momentum variable vanishes. Thus, the singular 
planes for each factor are: 

TTi = {p I Pi = 0}, (13.6) 

7r2 = {p\p2 = 0}, (13.7) 

^3 = {p\Pi-P2 = 0}, (13.8) 

i.e. vTi consists of all p such that pi = 0, etc. 

The set of all IR-subgraphs of G comprises three subgraphs 7^, i = 1,2,3, each 
consisting of one factor, and G itself. The two-factor subsets do not satisfy the com- 
pleteness condition and therefore are not IR-subgraphs. 



13.2 Fixing singularities corresponding to subgraphs 



Now we note that in the space of p without the point p = 0, the planes tTj do not 
intersect, so that a reasoning similar to that of subsects. |7.10| and |12.1| can be used. 
Indeed, consider a test function (pi{p) which is non-zero within the region O shown in 
Fig. 7b and therefore is zero in neighbourhoods of both 7C2 and ir^. Then 



where 



dpidp2G{p, K)(pi{p)= J dpigi{p, k)(pi{j)i,k) 



^i(Pi,/«) = / dp2G\gi{p,K)Lpi{p) = / dp2 g2ip, n) gsip, k) Lpi{p) 



[13.9) 



(13.10) 



The expansion of the expression |13.9| can be obtained (cf. subsect. |11.3| ; also note 
that we present the formulae in dimensionally regularized form) by Taylor-expanding 
<fi{pi, k) (which does not give rise to any singularities in pi) and using |7.40| to expand 
the propagator gi (cf. |1 1 . 1 1] ) . Denote: 



[Efi-l] (p, k) =J2 Ca{K,) SaiPl] 



:i3.iii 



with Cq,(k) given by [7.41| , where Sa and Va have been introduced in [7.36| , [7.37| (it is 
convenient to assume that 60 = 6 and V = 1). Then one can rewrite the obtained 
result as follows: 



G{p, k) ~ {Tgiip, k) g2{p, k) g^{p, k)) + {Tg2{p, g^ip, k)) [Egi] {p, k), (13.12) 

which is valid on the test functions (pi described above. 

Analogous expansions can be obtained for the cases of test functions ip2 and (fs that 
are in the same relation to the singular planes 7^2 and 713, respectively, as is to vri. 
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Consider a test function foip) which is identically zero in some neighbourhood of 
the point p = 0. It can always be represented as 



(13.13) 



where all ipi are as described above. Therefore, the three expansions that are valid on 
each Lfi can be glued together into an expansion that is valid on any ip^ : 



Gip, k) 



(T(7i(p, k) g^ip, k) g^ip, k)) + iTg2{p, k) g^{p, k)) [E(7i] (p, k) (13.14) 



+ (Tc/i(p, k) c/3(p, /«)) [E^2] (p, 1^) + (Tc/i(p, k) c/2(p, /«)) [E^a] (p, «;) 
= As'G{p, k). 
This agrees with |12.1| - |12.^ . 

Consider the expansion |13.14| to order o{n'^) : 

1 



(13.15) 
(13.16) 



As'^Gip, k) 



+ Cq[k) 
Recall that 



— K 



S{Pi 



Pipl{.pi-P2y 

1 1 

+ 



+ 



pipI (pi -P2y P1P2 {.P1-P2Y P1P2 {.P1-P2 



+ ^KP2) 



+ 



^2(^1-^2)^ pi {.pi-P2f pipl 



KPI-P2) 



Co(/«) = j d^q (g' + ft 



2\-l 



K 



2(D-4) 



Co(l). 



(13.17) 

(13.18) 
(13.19) 

(13.20) 



One sees that the O(k^) terms possess a logarithmic singularity at p = 0. Therefore, 



the r.h.s. of |13.17| is well-defined on test functions (y9(p) such that <^(0) = 0. Moreover, 
on such test functions the asymptotic character of the expansion is preserved which 
means (to order o{k^)) that: 



j d'^'p^ip) [Gip, k) - As',Gip, k)] = o{ 



(13.21) 



This should be compared with \l.3[ An explicit proof of |13.21| proceeds along the same 
lines as that of |7]^. A new element here is that one should consider cone regions in the 
integration space, as follows. (There is a limited resemblance between such cone regions 
and the Hepp sectors used in the theory of a-representation. However, the Hepp sectors 
correspond to a complete resolution of the recursive structure of singularities, while in 
our case only one-level descent — form the complete graph to its largest subgraphs — is 
performed.) 
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13.3 Approximation properties of As' 

Let 9i{p) be such that 6i{p) = 1 around each non-zero point of vti, 9i{p) = around 
non-zero points of it2 and tt^, and Oi{Xp) = Oi{p) for each X, p ^ 0. Consider the 



following contribution to 13.21: 



d'^pe,{pMp) [G{p, k) - As',G{p, «:)] . 



;i3.22) 



One splits the integration region into two parts: \pi\ > k, and |j9i| < k, (cf. Fig. 7c and 
subsect. [7.2|) . In the region < k, one rescales p ^ up then takes into account the 
scaling properties of G (which are not violated by the As-operation) and the fact that 
i^{p) ~ p near p = 0, and finds that 



|pi|<K 



d'^'pe.ip) ^{p) [G{p, k) - As',G{p, k)] = 0{k' 



fl3.23) 



In the region \pi | > it is convenient to reorder the integrand of |13.22| as follows: 

1 



Pl|>K 



d'''pe^{p)^{p){G{p,K) 



+ 



+ 



ViVl (Pi - P2 ) ^ P1P2 {.P1-P2Y P1P2 {V1-V2Y 



Pl|>K 



d'^peM^{p) X 



+ 



" 1 ■ 




_pI_ 





1 



2 I 9 1 ~^ ^ A 

PI + PI p{ 



1 



p\ + {pi - P2Y + ^2 
/ 1 1 



(13.24) 

(13.25) 
(13.26) 



pI + {pi - P2Y + ^2 pl{pi - P2Y 

1 1 

+ 



+ 



21" 








Pi. 





pI + N? {pi - P2Y + pUpi - P2y 
Consider e.g. the first term contributing to the r.h.s.: 



,P2{Pl~P2y P2{Pl-P2Y 

(13.27) 
(13.28) 



d^'>pe,{p)^{p) 



1pi|>k 



X 



1 



1 



Pi + K'^ Pi 



^ + K ^ - Co{k)8{Pi 



Pi 



P2 + {pi - P2Y + 



(13.29) 

(13.30) 
(13.31) 

■{p2 + n^){{pi-P2? + K^)' ^^^'^^^ 

The square-bracketed term on the r.h.s. of p.3.29| is bounded by const x /\pi\^. To 
estimate one should take into account that the integration over p2 in |13.32| runs over 



|pi|>K 



d°Piip{pi, k) 



2 I 9 2 ~^ ^ A 

Pi + K"^ Pi Pi 



where 



(^(Pi,k) = j d^p2- 
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the region \p2\ < k\pi\, and that ip can be represented as pi^pi + P2'^2, where both ipi 
are smooth functions Rescahng the integration variable as p2 — > \pi\-P2, one gets: 



I5 



iPil 



|D~4+1 



d 



P2\<k 



P2 



OliPl,P2)VliPl, 


Pi 


■P2) 


{pI + KVbiP)((pi -P2y 







(13.33) 



Since \pi \ ■ P2) can be bounded by const x \pi\, one sees that <^(pi, k) is bounded 

by const x |pi | . (It is not difficuh to understand that since the expansion |13.14| should 
be asymptotic only at Z) = 4, one may check estimates like |13.21| only a.t D = 4.) One 
arrives at the conclusion that the expression |13.29| is o(k^), as expected. The remaining 
terms in |13.24| are considered in the same way and, finally, one obtains |13.21, 



13.4 The last counterterm 



The result of the preceding subsection (the estimate |13.21| ) allows one to repeat the 
reasoning of subsect. |12.3| . 

One finally obtains: 

G{p,k) ~ As2G{p,k) +o{k'^), (13.34) 

where 

As2G(p, k) = As'^Gip, k) + ccMSip) (13.35) 

{5{p) = 5{pi)S{p2)) and 

cg,o(«:) = / d^''pG{p,K) = J d^p^d^p2GipuP2,f^)- (13.36) 
Examples of higher-order counterterms are: 

CG,2,o{^)d'6{p,)Sip2), (13.37) 

where 

cg,2,o(«) = J d'^Pi d'^P2Pl G{pi,p2, n), (13.38) 

and 

CGMd^Sipi)d^'Sip2), (13.39) 

where 

CG,o,2ifi) = J d'^Pi d^P2 {Pi ■ P2)G{pi,p2, k), (13.40) 
etc. Concerning diagrammatic interpretation of |13.36| , |13.38| and |13.40| , see Fig. 7d. 
This completes our discussion of the example. 

The essential point to be stressed here is that the proof involves nothing but the 
power counting. The same remains true in the most general case as well, which justifies 
the conclusions of subsect. |12.2| . For a more formal treatment of this point see |]20[ . 
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Note that in practice there is no real need to perform the detailed analysis as above. 
It is simply sufficient to enumerate singularities of the formal expansion, determine (by 
power counting) what counterterms are needed, and to find counterterms by formal 
integration of the ansatz for the ^s-operation with suitable polynomials to project out 
the coefficients. 
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Figure captions 
Fig.l. Kinematics of tliree-point functions. 

Fig. 2. (a) Kinematics of light-by-light scattering. It differs from the general 4-point 
function (b) by the restriction qi + qs = 0. 

Fig. 3. Three light-mass propagators connected to a common vertex form a complete 
IR-subgraph. 

Fig. 4. The property of being an IR-subgraph depends on restrictions imposed on 
heavy external momenta: in (a) the pair of vertical lines constitutes an IR subgraph, 
while in (b) it does not. 

Fig. 5. This graph is its own IR-subgraph if the contact terms localized at Q = 
are to be taken into account; it is not, if Q is just fixed at a non-zero value. 

Fig. 6. A simple example of the division of the proper variables of IR-subgraphs 
into "internal" components corresponding to loops and the " external" components. 

Fig. 7. (a) A graphical representation of the product |13.2| . (b) The test function ipi 
is non-zero within the region enclosed by the dashed line, (c) Geometry of integrations 
13.29| and |13.34 (d) A graphical interpretation of the expressions |13.36| , |13.38| and 



m 



13.40. 



